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Inclusive transverse momentum spectra of r| mesons in tiie range p/- ~ 2 - 12 GeV/c have been measured 
at mid-rapidity (|ri| < 0.35) by the PHENIX experiment at RHIC in p+p, d+Au and Au+Au collisions at ^JY^ 
= 200 GeV. The T] mesons are reconstructed through their T| ^ YY channel for the three colliding systems as 
well as through the r\ Jt^Tt+TT^ decay mode in p+p and d+Au collisions. The nuclear modification factor in 
d+Au collisions, RdAuipr) ~ 1-0 - 1.1, suggests at most only modest px broadening ("Cronin enhancement"). 
In central Au+Au reactions, the T] yields are significantly suppressed, with RauAu{pt) ~ 0.2. The ratio of r| to 
7t*' yields is approximately constant as a function of pj for the three colliding systems in agreement with the 
high-pr world average of Ry^/jfi ~ 0.5 in hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions for a 
wide range of center-of-mass energies {^/s^ ^ 3 - 1800 GeV) as well as, for high scaled momentum Xp, in 
e+e^ annihilations at \/s — 91.2 GeV. These results are consistent with a scenario where high-pj- T| production 
in nuclear collisions at RHIC is largely unaffected by initial-state effects, but where light-quark mesons 
are equally suppressed due to final-state interactions of the parent partons in the dense medium produced in 
Au+Au reactions. 
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I. INTRODUCTION 

Single-hadron production at large transverse momenta 
(Pt ^ 2 GeV/c) in high-energy hadronic and nuclear colli- 
sions results from the fragmentation of quarks and gluons 
issuing from parton-parton scatterings with large momentum 
transfer Since the cross-sections for such hard processes 
can be calculated perturbatively within Quantum Chromo- 
dynamics (QCD) flj, inclusive high-pj hadrons (as well as 
jets, real and virtual direct photons, and heavy quarks) have 
long been considered sensitive, well calibrated probes of the 
small-distance QCD processes. The study of inclusive hadron 
production at large pr in proton-proton interactions provides 
valuable information about perturbative QCD (pQCD), parton 
distribution functions (PDF) in the proton, and fragmentation 
functions (FF) of the partons |2]. Furthermore, the use 
of polarized beams {p + p) allows one to investigate the 
spin structure of the proton |3]. High-energy collisions of 
protons or deuterons on nuclear targets (p,d+A) also provide 
interesting insights on initial- and final-state QCD effects 
such as modifications of the nuclear PDFs [4, 5] and parton 
rescattering in the cold nuclear medium Both effects are 
sensitive to physics such as parton structure and evolution 
at small values of fractional momentum x in the hadronic 
wave functions JED, and the dynamics of hadronization in 
cold nuclei JH Lastly, high-pj- hadron production in 
nucleus-nucleus (A+A) reactions is a sensitive probe of the 
properties of the strongly interacting QCD matter produced 
in the collision. Indeed, since perturbative processes happen 
at time scales T w 1 /f r ^ 0. 1 fm/c, the hard-scattered partons 
traverse and are potentially modified by the bulk matter 
formed shortly after the collision. In this context, the suppres- 
sion of leading hadrons has been postulated ITdh as a signal 
of "jet quenching" in a Quark-Gluon-Plasma (QGP) due to 
medium-induced energy-loss of the parent parton 1 11 , 1 ^ [Tsl] . 

All the aforementioned research topics have been addressed 
in detail by the rich physics program carried out at the Rela- 
tivistic Heavy-Ion Collider (RHIC) at Brookhaven National 
Laboratory (BNL) during its first six years of operation (2000 
- 2006). For example, the study of inclusive high-pj- neu- 
tral pion production at mid-rapidity in p+p 13, p + p ifisll . 
d+Au [16, 17] and Au+Au [18, 19,i3,lMl collisions at ^/i^ 
= 200 GeV, has provided valuable information respectively 
on: 

(i) the gluon-to-pion FF lHHH, 

(ii) the gluon contribution to the proton spin ll23il24ll25ll . 

(iii) initial-state effects in cold nuclear matter such as shad- 
owing of the nuclear PDFs IH HUll IH, Cronin 
broadening [HI HI HI, and gluon saturation IH [M 
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(iv) the thermodynamical and microscopic properties of 
hot and dense QCD matter liS^]. such as the initial 
gluon rapidity density dN^ /dy 0371] and the transport 
coefficient {q) ifssl l39ll of the produced medium; and 
the mechanism of hadronization in a dense parton 
medium ll40ll . 



In this paper, we extend previous PHENIX analyses of 
high-pr hadron production in p+p ifTil [Tsll . d+Au |16] and 
Au+Au 1 19, 201 111 SH colHsions at = 200 GeV, to in- 

clude an additional identified particle, the r| meson, measured 
in the range 797- = 2 - 12 GeV/c. The spectra presented here 
are the hardest (i.e. have the highest pr) ever measured for the 
r| meson' in p+p, p+A and A+A collisions. The high pr reach 
of the r| helps to characterize the mechanisms of truly pertur- 
bative parton-parton scatterings and parton fragmentation in 
different QCD environments (p+p, d+A and Au+Au). The rj 
data from p+p collisions are presented here as a baseline for 
medium effects in d+Au and Au+Au. Once a parametriza- 
tion of the r| FF in e+e^ is performed (see Section IIVE2I ). 
the observed p+p cross section will additionally allow a test 
of pQCD predictions. Such a FF parametrization would be 
useful in particular in the light of upcoming high-pj- r| asym- 
metry data obtained with polarized beams of relevance for the 
proton spin program [44]. 

For d+Au and Au+Au reactions, we present the single r| 
spectra, the r| nuclear modification factors, and the rj-to-Tt" 
ratio measured as a function of pr in different centralities. 
Within uncertainties, the d+Au spectra for all centralities are 
consistent with the p+p yields scaled by the corresponding 
number of incoherent nucleon-nucleon (A^A^) collisions. The 
maximum amount of pj broadening seen in the r| data is 
10%, RdAtiiPr) ~ 1-0 - 1.1. Such a result confirms the 
limited influence of cold nuclear matter effects, such as 
shadowing, Cronin broadening or recombination, on high-p^ 
meson production at mid-rapidity at RHIC [|16l llTll . On the 
other hand, the factor of ^5 deficit of inclusive r| yields 
observed above pjK, A GeV/c in central Au+Au compared to 
binary-scaled p+p collisions, RauAu{pt) ~ 0.2, is the same as 
that found for YAgh-pj 1, 19, 21] and for inclusive charged 
hadrons ll4lll45ll . Such a common suppression pattern for Tt", 
r| and /1* is expected if the energy loss takes place at the par- 
ton level in the dense medium formed in the reaction prior to 
its fragmentation into a given hadron in the vacuum. Indeed, 
in this case the high-pj- deficit will just depend on the energy 
lost by the parent light-quark or gluon (i.e. on the initial 
density of scatterers in the produced medium) and not on the 
nature of the final leading hadron whose production will be 
described by the same universal probabilities (fragmentation 
functions) which govern vacuum hadron production in more 
elementary systems. Such an interpretation is supported 



' Before this measurement, only the ISR AFS Collaboration p+p^ r|+X 
measurement for = 3 - 1 1 GeV/c a\ \/s = 62.4 GeV [42] and the single 
CDF ri/rc" point measured at pj = 12 GeV/c in p + p collisions at. \/s = 
1800 GeV 1 43] had comparable maximum pj values. 



by the fact that the ratio of /jj-differential cross-sections 
of r| mesons with respect to ii^ in Au+Au, d+Au and p+p 
collisions is approximately constant, R^it^o ~ 0.40 - 0.48, 
which is consistent with the world average measured: (i) in 
hadron-hadron, hadron-nucleus and nucleus-nucleus colli- 
sions above pr ^ ^ GeV/c; as well as (ii) in electron-positron 
annihilations at the Z pole =91.2 GeV) for energetic 
r| and 71'' with scaled momenta Xp = Phadmn/ Pbeam ^ 0.35. 
Comparison of our data to a world compilation of rj/Tt" ratios 
is done in the last section of the paper. 

In addition to their interest as a signal in their own right, 
reliable knowledge of the production of r| mesons in p+p, 
d+Au and Au+Au reactions is also required in order to deter- 
mine and statistically remove the background of secondary 
and y for other measurements such as single electrons 
(from heavy-quark decays) ll46l I47I1 . dielectrons ll47ll and 
direct photons [48, "49"]. Indeed, r| mesons constitute the 
second most important source after the Tt" of decay electrons 
(Dalitz and conversion) and photons contributing to these 
backgrounds. 

The paper is organized as follows. Section |ll] presents a 
description of the experimental setup and detector systems 
used in this work. Section |III] provides an explanation of 
the analysis methods employed to obtain the r| data. Sec- 
tion |IV] presents and compares the r| results (pr spectra, nu- 
clear modification factors, and ratios) measured in p+p, 
d+Au and Au+Au at ^/s^ = 200 GeV, and discusses the rel- 
ative role of cold nuclear (d+Au) and hot and dense medium 
(Au+Au) effects on high-pr meson production. In particular, 
Section lTV El discusses the measured rj-to-Ti" ratios in the con- 
text of different phenomenological models of high-pj- hadron 
production as well as in comparison to other experimental re- 
sults measured in high-energy particle collisions at different 
center-of-mass energies. A less detailed presentation of a sub- 
set of these r| results has already been published in ll20tl . 

II. EXPERIMENTAL SETUP 

The PHENIX experiment at the RHIC facility lH at BNL 
is specifically designed to measure hard QCD probes such 
as high-p7 hadrons, direct photons, leptons, and heavy flavor 
production. PHENIX achieves good mass and particle identi- 
fication (PID) resolutions as well as small granularity by com- 
bining 13 detector subsystems (^350,000 channels) divided 
into: (i) two central arm spectrometers for electron, photon 
and hadron measurements at mid-rapidity (|ri| < 0.35, A(|) = 
71); (ii) two forward-backward (|r| I = 1.2 - 2.2, A(|) = 2%) spec- 
trometers for muon detection; and (iii) two global (inner) de- 
tectors for trigger and centrality selection. A detailed descrip- 
tion of the complete detector can be found elsewhere Iplll . 
The data presented in this paper were obtained during the 
Run-2 (2001-2002, Au+Au) and Run-3 (2003, d+Au, p+p) 
operations at RHIC. The layout of the PHENIX detector dur- 
ing these run periods is shown in Fig.[T] The primary detectors 
used to obtain the present results are the PHENIX central arm 
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spectrometers, particularly the electromagnetic calorimeters 
(EMCal) 1 52] and the charged particle tracking devices (the 
Drift Chamber (DC) lH and Pad Chambers ^C) Q). 
In addition, the Beam-Beam Counters (BBC) I5y\ and the 
Zero-Degree Calorimeters (ZDC) |56] are used for triggering, 
event characterization and (AuH-Au and dn-Au) centrality 
determination. 




FIG. 1 : The PHENIX experimental setup during Run-2 and Run-3 
at RHIC. The detectors used in the present analysis are the 8 EMCal 
(PbSc, PbGl) sectors for photon detection (rj yy), the drift cham- 
ber (DC) and two layers of multiwire proportional chambers with 
pad readout (PC) for charged pion detection (r| 7t"7l+7t^); as well 
as the two Beam-Beam Counters (BBC) and the two Zero-Degree 
Calorimeters (ZDC) for global event characterization. 



A. Electromagnetic Calorimeter (EMCal) 

The r| mesons are detected in PHENIX via their yy 
(branching ratio BR = 39.43%) and Tz\+%- (BR = 22.6%) 
decays 115711 . Photons from the direct yy channel as well as 
from the secondary (daughter) Tt** decays are measured in the 
PHENIX EMCal, which has a pseudo-rapidity acceptance 
of —0.35 < T] < 0.35 and covers K radians in azimuth. The 
electromagnetic calorimeter is divided into eight sectors with 
two distinct detection technologies (see Fig. [T]). A lead- 
scintillator calorimeter (PbSc) consists of 15552 individual 
lead-scintillator sandwich modules (5.54 cm x 5.54 cm x 



37.5 cm, ISXq), grouped in six sectors located at a radial 
distance of 5.1 m from the beam line, covering a total solid 
angle of Ar| w 0.7 and A(|) « 37t/4. A lead-glass Cerenkov 
calorimeter (PbGl) comprising two sectors, with a total of 
9216 modules (4 cm x 4 cm x 40 cm, IAAXq), is located at 
a radial distance of ^5 m from the beam pipe and covers a 
total solid angle at mid-rapidity of Ar| sa 0.7 and A(|) = 7t/4. 
The corresponding Ar| x A(|) acceptance of a single tower at 
ri = is 0.011 X 0.011 and 0.0075 x 0.0075 for the PbSc 
and PbGl calorimeters, respectively. The chosen transverse 
size of the towers is not much larger than their corresponding 
Moliere radius (pM = 3.0 cm and 3.7 cm for PbSc and PbGl, 
respectively) so that most of the electromagnetic showers 
extend over several modules resulting in an improved position 
resolution based on a "center of gravity" reconstruction of the 
impact point of the photon clusters. 

The energy calibration of the PbSc modules was obtained 
from the original beam-test values and redundantly confirmed 
with (i) the position of the reconstructed Tt** mass peak, (ii) 
the energy deposit from minimum-ionizing charged particles 
traversing the calorimeter, as well as with (iii) the expected 
EphScI Ptrucking ~1 value measured for electrons and positrons 
identified in the Ring-Imaging Cerenkov (RICH) detector and 
whose momentum was measured in the tracking detectors. In 
the PbGl modules, the reference energy calibration from the 
original beam-test values is corrected with the time-dependent 
gain corrections obtained with a light-emission-diode (LED) 
system for the lead-glass calorimeter. The LEDs emit light 
with known intensity, so gain fluctuations can be detected. 
The final PbGl calibration is obtained by comparing the mea- 
sured Tt'' peak position to its nominal value. 



B. Central Arm Tracking 

Charged pions are measured with the PHENIX central 
tracking system combining information from the drift and 
pad chambers. The momenta of the 7t* are measured at a 
radius of 2.0 m from the event vertex by the Drift Chamber 
(DC). The DC, located outside the field of PHENIX central 
magnets, uses several layers of wires to reconstruct the 
angle of the track, which is inversely proportional to its 
momentum. The DC momentum resolution is determined 
to be Q.1®\.\%pt (GeV/c). The polar angle of the track 
is measured by Pad Chamber 1 (PCI), a multi-wire pro- 
portional chamber located just beyond the DC. The last 
pad-chamber layer, PC3, at a radius of 5.0 m and directly in 
front of the EMCal, is used in this analysis for two purposes: 
to confirm the track by matching to a PC3 hit, as well as 
to veto an EMCal cluster produced by a charged particle track. 

The DC momentum scale is checked by the reconstruction 
of the correct mass of (i) 7t*, K^, p,p identified with the time- 
of-flight system |58], and (ii) CO, (|), J/^ mesons decaying into 
the e+e" channel identified with the RICH and EMCal. The 
momentum scale is thus known with an accuracy better than 
0.2%. Since at low pj, the momentum resolution is better 
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TABLE I: Events sampled and integrated luminosity (after vertex cuts) in the T] analyses for p+p, d+Au and Au+Au collisions. The equivalent 
p+p luminosities in d+Au (Au+Au) have been obtained normalizing their corresponding luminosities by 2A (A^) factors as expected for hard 
cross-section scaling. 



collision 


Events sampled 


Total integrated luminosities 


BBC attributes 


system 


MB (LVLl) trigger 


High-pj- (Gamma) trigger 


absolute 


equivalent p+p 


cross section efficiency 


p+p 


25.2- 10'' 


49.3-10** 


216 nb"^ 


216 nb"' 


23.0 mb ± 9.7% (55±5)% 


d+Au 


58.3 10^ 


29.2- 10^ 


1.5 nb^i 


590nb-l 


1.99 b ±5.2% (88±4)% 


Au+Au 


34- 10^ 


30 - lO*" 


9ijb-^ 


230nb-i 


6.315 b ±8.4% (92±3)% 



when measured with the tracking system than that using the 
energy measured via calorimetry, and given that the momen- 
tum range of the three r| decay products has a relatively low 
Pt, the uncertainties in the tracking system calibration are less 
important in the Tt^Ti+Tt^ measurement than in the yy decay 
channel. As a result, the tracking devices provide better accu- 
racy for the Tj mass reconstruction than the EMCal. 



C. Global Detectors 

Triggering and global event characterization is achieved 
using the Beam-Beam Counters (BBC) and the Zero-Degree 
Calorimeters (ZDC). The two BBC are placed around the 
beam pipe 1.44 m in each direction from the nominal in- 
teraction point. Each BBC consists of 64 hexagonal quartz 
Cerenkov radiators closely packed in an approximately 
azimuthally symmetric configuration. The BBC are used to 
count the charged particle multiplicity in the pseudo-rapidity 
range 3.0 < |r|| < 3.9, to provide the start time for time-of- 
flight measurement, and to give the collision vertex position 
along the interaction diamond with a typical resolution of 0.6 
(2) cm in Au+Au (p+p) collisions fSSll . In d+Au collisions, 
the centrality of the collision is determined by measurin g th e 
charge deposited in the BBC in the Au beam direction 115911 : 
whereas in Au+Au reactions, the correlation between the 
BBC charge sum and the ZDC total energy is used for 
centrality determination |60] (see next Section). The ZDC are 
small hadronic calorimeters that measure the energy carried 
by spectator neutrons at very forward angles. They are placed 
18 m up- and downstream of the interaction point along 
the beam line (|9| < 2 mrad). Each ZDC consists of three 
modules with a depth of 2 hadronic interaction lengths read 
out by a single photo-multiplier tube (PMT). Both time and 
amplitude are digitized for each PMT along with the analog 
sum of the three PMT signals for each ZDC. 



III. DATA ANALYSIS 

In this section, we describe the event selection criteria, the 
reaction centrality determination in d+Au and Au+Au colli- 
sions, the r| identification and reconstruction procedures in the 
Tj ^ YY and r| — > ti^ti+ti^ channels and the various correc- 
tions (geometrical acceptance, reconstruction efficiency, trig- 
ger, absolute cross-section normalization) applied to the raw 



data. The systematic uncertainties of the measurements are 
discussed at the end of the section. 



A. Event Selection 

The data presented in this paper were collected under two 
general trigger conditions. The first sample, consisting of 
minimum-bias (MB) events with vertex position along the 
beam axis |z| < 30 cm, was conditioned on a local-level- 1 
(LVLl) trigger generated by coincidences between the two 
BBC (in the case of p+p and d+Au) or by coincidences be- 
tween the BBC and ZDC detectors (in the case of Au+Au). 
The MB trigger cross sections measured by the BBC in p+p 
and d+Au collisions are respectively 23.0 mb ± 9.7% mb 
and 1.99 b ± 5.2% |61], whereas the Run-2 Au+Au mini- 
mum bias trigger has some inefficiency for the most periph- 
eral interactions and records only 92.2^3 q% of OauAu 119(1 . In 
other words, the LVLl triggers accept respectively (55±5)%, 
(88±4)%, and (92±3)% of the total inelastic cross sections: 
a'l^'p' = 42 ± 3 mb, a'X^ = 2260 ± 100 mb, and a^L = 6850 
± 540 mb. A second "photon-triggered" sample, requiring 
electromagnetic showers above a given threshold in the EM- 
Cal (with or without the MB BBC requirement), has been used 
to extend the r| measurements to higher pr- The details of this 
level-2 (LVL2) software trigger are described in li2lll . The to- 
tal number of events collected in the MB and photon-triggered 
samples (after vertex cuts) as well as the integrated luminosi- 
ties for each collision system are listed in Table U 



B. Centrality Determination (d+Au, Au+Au) 

The events in d+Au collisions are classified into four dif- 
ferent centrality classes given in percentiles of the total cross- 
section: 0-20%, 20-40%, 40-60% and 60-88%, with the lat- 
ter being the most peripheral. The reaction centrality is re- 
lated to the number of hits in the south Beam-Beam Counter 
(BBCS), which is proportional to the number of participant 
nucleons in the gold nucleus ll59ll . The distribution of the nor- 
malized charge in the BBCS and the classification into differ- 
ent centrality classes is shown in Fig. |2l In order to obtain rea- 
sonably large statistics in each Au+Au centrality class, three 
centralities are used in the current Au+Au analysis: 0-20% 
(central), 20-60% (semicentral) and 60-92% (peripheral), de- 
termined by cuts in the correlated distribution of the charge 
detected in the BBC and the energy measured in the ZDC Ii60fl . 



7 



Ii3 

o to 
c m 

o 



9000 
8000 
7000 
6000 
5000 
4000 
3000 
2000 
1000 





FIG. 2: Distribution of tiie normalized ciiarge in the south Beam- 
Beam Counter (BBCS) in d+Au collisions at ^/s^ = 200 GeV. The 
normalization is done such that the normalized charge corresponds 
to the number of hits. 



A Glauber Monte Carlo model combined with a simulation 
of the BBC (plus ZDC) response allows determination of the 
mean value of the associated nuclear overlap function {Tja) 
({Taa)) for each d+Au (Au+Au) centrality bin. Table lists 
the mean value of the nuclear overlap function for different 
centralities in both systems. 



TABLE II: Values of the average nuclear overlap function {T^a} ^nd 
(Taa) for the different centralities considered in d+Au and Au+Au 
reactions respectively. 



centrality bin 



(r,M) (mb-') 



(Taa) (mb-') 



min bias 
0-20% 
20-40% 
20-60% 
40-60% 
60-88% 
60-92% 



0.20 ± O.OI 
0.36 ± 0.02 
0.25 ± 0.017 

0.17 ± 0.014 
0.073 ± 0.007 



6.14 ±0.45 
18.5 ± 1.3 

4.6 ± 0.4 



0.3 ± 0.1 



C. r| YY Reconstruction 

The main mode of r| -meson reconstruction in PHENIX is 
via the electromagnetic channel r| yy. PHENIX has pub- 
lished the results of a number of 7t*^ yymeasurements in 
the EMCal for different colliding systems 11141 
[l9[ I201 12? ] . The technique for identifying the photons and 
reconstructing the 7i'' yields as a function of pr and central- 
ity is now well established and is exactly the same one used 
here to obtain the corresponding r| yields. Although the recon- 
struction methods are identical, the p+p and d+Au analyses do 
not suffer from the large particle multiplicity background that 
the Au+Au r| reconstruction faces, and there are a few differ- 
ences between Au+Au and the other studies. In the Au+Au 
case, the applied analysis cuts (photon identification, invari- 
ant mass reconstruction, and other cuts) are tighter than in the 



p+p and d+Au cases. Additionally, in order to deal with clus- 
ter overlap effects appropriately, the Au+Au analysis uses a 
full GEANT [63] Monte Carlo (MC) simulation, in which the 
simulated single r| are embedded into real events, rather than 
using a tuned fast MC simulation without embedding. These 
differences are explained in separate subsections below. 



1. Photon reconstruction in EMCal 

Electromagnetic clusters are reconstructed in the EMCal 
sectors by finding contiguous calorimeter towers with pulse 
heights above the ADC pedestal value. The energy of each 
EMCal cluster is corrected for angular dependence and 
non-linearity based on test-beam results and simulation. The 
linearity correction for the two detector types is different, 
with the PbGl having a stronger dependence on the energy. 
The correction factors for a photon with a detected energy 
of 1 GeV (10 GeV) are 1.00 (0.95) for the PbSc and 1.05 
(0.975) for the PbGl, respectively. The PbGl calorimeter 
also shows a stronger variation of the measured photon 
energy with the angle of incidence on the detector surface: 
at 20° the measured energy is reduced by 5% compared to 
perpendicular incidence (0°), while in the PbSc this reduction 
is only of the order of 2%. 

Since we are interested in high- pr T| production, only 
EMCal clusters with energy above 1 GeV are selected for 
further analysis. In addition, fiducial cuts excluding the 
edges of each EMCal sector as well as an area of 3 x 3 
towers around the towers that have been determined to be 
hot or dead, were applied in order to exclude clusters with 
incorrectly reconstructed energies. Among the clusters 
passing the cuts, photon candidates are identified by applying 
standard particle identification (PID) cuts based on time- 
of-flight (TOE) and shower profile. Both cuts are applied 
to reject slower and broader showers which are mostly of 
hadronic origin. For the PbSc we require the measured cluster 
TOE to be fciust < L/c± 1.2 ns where L is the straight-line 
path from the collision vertex to the reconstructed cluster 
centroid. For the PbGl we require reconstructed clusters to 
have times, fdust < L/c±2 ns; the difference is due to the 
difference in intrinsic timing resolution of the two calorimeter 
technologies. Shower profile cuts are based on rejecting those 
clusters whose energy deposition among the modules, and 
in particular in the most central tower of the cluster, are not 
consistent, within a given x^-test limit, with the shower shape 
expected for electromagnetic showers as parametrized from 
test-beam data 115 2[]. 



In the most central Au+Au events, the EMCal typically de- 
tects >300 clusters, corresponding to a detector occupancy 
of ^10% in terms of hit towers, resulting in a non-negligible 
probability of particles making clusters which overlap. In or- 
der to minimize the effects of cluster overlaps due to high 
multiplicity, two methods are used to determine the cluster 
energy. First the energy of each cluster in the PbSc calorime- 
ter is determined from the sum of all contiguous towers with 
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deposited energy above a given threshold (Etowei- =15 MeV, 
typically). Alternatively, an extrapolation (using a standard 
electromagnetic shower profile for an event with zero back- 
ground) from the measured core energy (ecore) in the four 
central towers to the full cluster energy is used. For the latter 
case, the ecore energy was computed from the experimentally 
measured center of gravity, central shower energy, and impact 
angle in the calorimeter using a parametrized shower profile 
function obtained from electromagnetic showers measured in 
the beam tests. Such an ecore energy represents an estimate of 
the true energy of a photon impinging on the PbSc unbiased 
by background contributions from other particles produced in 
the same event and depositing energy in the neighborhood of 
a given cluster. The use of ecore instead of the total cluster 
energy for photon reconstruction helped considerably to min- 
imize the effects of cluster overlaps in central AuH-Au colli- 
sions. 



2. Raw T| yield extraction (p+p and d+Aii) 

The r| yields are obtained by an invariant mass anal- 
ysis of photon candidate pairs having asymmetries 
a = \Ey^ —Ey^\/ {Ey^ + Ey^) < 0.7. Thc cut on the asymmetry 
a reduces the background since high-pj- combinatorial pairs 
are strongly peaked near a = 1 due to the steeply falling 
spectrum of single photon candidates. The resulting invariant 
mass spectra obtained for proton-proton and deuteron-gold 
collisions are shown in Figs. [3] and |4] respectively for two 
typical pt bins. A peak is seen at about 550 MeV/c^, the 
expected mass of the r| meson. The measured peak position is 
modified by detector effects which lead to energy smearing. 
The combined effects of the energy resolution of the detector, 
the steeply falling single photon spectrum, and the finite size 
of the energy bins lead to a smearing of the measured photon 
energies which widens the r| signal. As a result the average 
peak position in the invariant mass spectra is about 9 MeV/c^ 
larger than the nominal mass of the r| meson, an effect which 
is well reproduced by the simulation. 

The combinatorial background below the peak signal is 
estimated with the event mixing method in which clusters 
from different events with similar event vertex (and centrality 
class in dn-Au) are combined to produce a "background" 
invariant mass distribution. This background is normalized 
to the real invariant mass spectrum and then subtracted from 
the invariant mass spectrum of the real events. To estimate 
the normalization of the background, the distribution of the 
real events is first divided by the mixed event distribution. 
This ratio is shown in the upper panel of Figs. [3] and |4] The 
normalization function is estimated by a fit in the region 
outside the peak. The spectrum, fitted to a second-degree 
polynomial, is shown in the region denoted by the vertical 
lines in the upper and the middle panels of Figs. [3] and ID 
The final real event mass distribution after the background 
subtraction is shown in the lower panel of Figs.[3]and|4l 

The interval over which the background is adjusted is lim- 
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FIG. 3: (color online) Invariant mass distribution of photon candi- 
date pairs measured in p+p collisions for the default PID cuts with 
pair transverse momenta 4.0 GeV/c < pj < 4.5 GeV/c. Top: Ra- 
tio of real and mixed event distributions, and background fits. The 
red fit is used for the background parametrization; the green fit for 
estimating the systematic uncertainty. Middle: Real invariant mass 
spectrum and scaled background. Bottom: Final distribution with the 
scaled background subtracted from the real event distribution (black 
entries); the green entries result from the background fit for estimat- 
ing the systematic error. Additionally, the peak is fitted with a Gaus- 
sian to get its mean and sigma. 



ited by two considerations: the expected r| peak position m 
and width o. Both were estimated in a first analysis of the 
spectra and set to m = 556 MeV /c^ and o 32 MeV /c^. The 
background interval includes the region between m — 7 .5o and 
m + 8.5o (320 MeV/c^ and 830 MeV/c^) excluding the peak 
region m ± 3a (460 MeV/c^ < m,„„ < 650 MeV/c^). For 
higher transverse momenta, the background almost vanishes 
and thus the estimation of the normalization by a fit leads to 
large errors. Hence an alternative method is used where the fit 
function is replaced by the ratio of the number of photon pairs 
in the normalization region in the real and the mixed event 
distributions. 

Finally, the total number of r| in a given pj bin is obtained 
by integration of the invariant mass distribution within 3o 
around the r] peak position. The statistical error of the peak 
extraction is estimated as done for the and described 
in |21]. The uncertainty of the background parametrization is 
estimated by calculating the error on the ratio of the integrals 
of the real and the mixed event distributions in the region of 
the background fit. Above pj = 10 GeV/c in pH-p collisions, 
the mixed event background does not work as expected as 
there are some entries in the mixed event background but 
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not in the r| region. For these cases, the background was 
estimated by integrating the real event distribution outside the 
peak (in the fit region mentioned above) and scahng this to 
the r| integration region. This background is then subtracted 
from the real event distribution to get the number of r|. The 
error in this case is estimated with o|,-^, = S + 2B, S being 
the r| signal and B being the background. The integral of the 
invariant mass distribution after the background subtraction is 
calculated in a reduced interval m — 2a; m + 2a (492 MeV/c^ 
- 620 MeV/c2). 
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FIG. 4: (color online) Invariant mass distribution of photon candi- 
date pairs measured in minimum bias d+Au collisions for the default 
PID cuts with pair transverse momenta 3.5 GeV/c < pr < 4.0 GeV/c. 
Top: Ratio of real and mixed event distribution, and background fits 
(the red fit is used for the background parametrization, the green fit 
for estimating the systematic uncertainty). Middle: Real invariant 
mass spectrum and scaled background. Bottom: Final distribution 
with the scaled background subtracted from the real event distribu- 
tion (black entries). The green entries result from the background fit 
for estimating the systematic error. Additionally, the peak is fitted 
with a Gaussian to get its mean and sigma. 



3. Raw T] yield extraction (Au+Au) 

The r| yields for AuH-Au are determined by calculating 
the invariant mass of photon pairs with asymmetries a = 
|£yi —E'^\/{Ey\ + E-fl) < 0.5, a value tighter than that used 
for the pH-p and dn-Au cases in order to reduce the larger un- 
corrected background in AuH-Au collisions, and binned in 
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FIG. 5: Top: Invariant mass distribution of pairs of photon candidates 
measured in minimum bias Au-l-Au with pair momenta pj = 3.5 - 
4.5 GeV/c around the r| mass fitted to a Gaussian plus exponential. 
Bottom: Final r[ signal after mixed-event (and residual) background 
subtraction. 



Pt- The r| yield in each pr bin is determined by integrat- 
ing the background-subtracted yy invariant mass distribution 
around the r| peak. The combinatorial background is obtained 
by combining uncorrelated photon pairs from different events 
with similar centrality and vertex, and normalizing the distri- 
bution in a region below (m,„v = 400 - 450 Me V/c^) and above 
(m,„„ = 750 - 1000 MeV/c2) the r] mass peak (Fig. |5] top). 
After the mixed background subtraction, the resulting distri- 
bution is fitted to a Gaussian plus an exponential (or linear, 
see below) function to account for the residual background - 
more important at low pr - not completely removed by the 
event-mixing technique. The bottom plot in Fig.|5]depicts the 
r| signal after mixed (and residual) background subtraction. 
To estimate the uncertainty in the subtraction of the residual 
background, different pair asymmetries and an alternative lin- 
ear parametrization were used (see Section lTlIFI ). The signal- 
to-background (S/B) ratio in peripheral (central) AuH-Au colli- 
sions is approximately 1.3 (1.5) and 0.05 (0.002) for the high- 
est and lowest pj, respectively. The signal-to-background ra- 
tio is comparable for central and peripheral collisions at the 
highest Pt because the spectrum in the central data extends to 
higher pj than that in the peripheral. 

The scaled mixed-event distribution is subtracted from the 
real-event distribution to produce a statistical measure of the 
true r| yield. The result of such a subtraction procedure is 
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shown in the bottom plot of Fig. |5] The raw r| yield is ob- 
tained by integrating the subtracted invariant mass distribu- 
tion in a range determined by the mean and the width of the r| 
peak and given by minv G [wri — 20n,'WTi +2ar|]. Th analysis 
described above is applied in bins of Apr = 1 GeV/c for pr 
= 2-4 GeV/c and Apr = 2 GeV/c above. We cease attempt- 
ing to extract r| yields at high pr when the number of pairs 
within the selected (background-subtracted) r] mass window 
falls below 4 counts. 



4. Acceptance correction 

The geometric acceptance is evaluated using a fast Monte- 
Carlo (fastMC) program based on routines from the JETSET 
library |64] to simulate the r| ^ yy decays and determine the 
geometric acceptance of the calorimeter for the decaying pho- 
tons. The acceptance correction accounts for the fraction of 
produced r| mesons whose decay photons will not actually hit 
the detector due to the finite solid angle covered by the de- 
tector A decay photon will be accepted by the EMCal in the 
fastMC when it hits the active surface of the detector covering 
the pseudorapidity range —0.35 < r| < 0.35 (computed using 
a realistic distribution of event vertices within |z| < 30 cm) 
and 2 x 90° in azimuth. The acceptance shows a strong depen- 
dence on the transverse momentum because the opening angle 
of the decay photons decreases with increasing pj- Thus, the 
probability that both decay photons hit the detector decreases 
for small values of pr- The acceptance for pH-p and dn-Au col- 
lisions is shown in Fig. |6] The acceptance is influenced by 
the geometry of the whole detector as well as by dead and hot 
modules in the pH-p and dn-Au cases (for AuH-Au the efficiency 
losses due to dead and hot modules are computed from the full 
GEANT3 simulation plus "embedding" and are accounted for 
in the efficiency loss correction). Due to a different number of 
masked out modules, the acceptance is not exactly the same 
in pH-p and dn-Au collisions. 



5. Efficiency corrections of the raw T) yields (p+p, d+Au) 

a. Reconstruction efficiency correction The reconstruc- 
tion efficiency takes into account that the measured r| spec- 
trum in the detector is different from the real physical spec- 
trum, i.e. the reconstruction efficiency is defined as the ratio 
between the output and the input r| spectra: 



dNr^/dprl 



out pill 



dNr^/dpT\i,ipui 



(1) 



as obtained using the fastMC, which parametrizes all the de- 
tector effects on the input spectrum (EMCal energy and po- 
sition resolution, efficiency losses due to y identification cuts 
and yy reconstruction procedure, etc.). A realistic input r| 
spectrum dN /dprlmput is used as an initial spectrum for the 
efficiency calculation and an iterative procedure is performed, 
in which the corrected output spectrum is used as the input 
spectrum of the next iteration. To simulate detector effects. 
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FIG. 6: Geometric acceptance (including dead channels) for the T| 
meson as a function of px measured in the EMCal in both d+Au 
(dashed curve) and p+p (solid curve) collisions at PHENIX in Run- 
3. 



the smeared energies and hit positions of the decay photons 
are parametrized in the fastMC. The energy smearing has a 
constant and an energy-dependent term and follows the func- 
tional form 



Oe/E- 



)B. 



(2) 



v/£/GeV 

The parameters for Eq. (|2]l are given in Table |III] for the 
different collision systems and the two EMCal detector types. 
The initial values have been taken from the detector response 
obtained in the beam tests ll52ll and re-tuned for real run con- 
ditions in previous iz^ analyses lfl4l ITsll . During the cross- 
checks between simulated and real data it was found that the 
energy scale of the EMCal was slightly shifted compared to 
the parametrized results. Since the energy scale is estimated 
experimentally by fitting the location of the Tt*^ peak, and the 
position of the peak is also affected by reconstructed sec- 
ondary neutral pions from K^^ decays which themselves decay 
off vertex, an additional correction is applied in the fastMC 
shifting the energy scale by 0.7%. After this correction, the 
position and the width of the simulated r] peaks are confirmed 
to be consistent with the position and the width measured in 
the data for all pj bins. 

TABLE III: Parameters for energy smearing, Eq. l|2j, as used in the 
fastMC for the different EMCal detector types and the different col- 
lision systems. 



collision 


detector 


energy-dependent 


constant 


system 


type 


term (A) 


term (B) 


p+p 


PbGl 


0.085 


0.059 


p+p 


PbSc 


0.082 


0.050 


d+Au 


PbGl 


0.085 


0.059 


d+Au 


PbSc 


0.082 


0.050 



The efficiency correction also takes into account the 
different cuts used for particle identification. The simulation 
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must consider the loss of photons and thus of r| due to the 
applied shower shape (or "dispersion") and energy threshold 
cuts. The effect of the dispersion cut is estimated by a 
comparison of uncorrected spectra without a PID cut with the 
spectra obtained with the different PID cuts. The spectra are 
obtained with a sharp asymmetry (a) cut and as a function of 
(£1 +E2)/2. The resulting loss of r| is translated into a photon 
loss probability, which is then used in the simulation. The 
energy cut is reproduced by rejecting photon hits according 
to an energy-dependent survival probability in the simulation. 
Finally, the simulation reconstructs the invariant mass and 
the transverse momentum of the r| from the reconstructed 
(smeared) information. Only particles inside the interval used 
for the integration of the real peak are accepted. The overall 
r| efficiency losses obtained by this method are of the order 
of Gri^YY = 76% ± 3% (dominated by the asymmetry cut and 
the invariant mass yield extraction procedure) and are flat 
within 1-2% in the whole pr range measured for both (p+p 
and d+Au) colliding systems. 

b. Photon conversion correction Some of the produced 
T] are not reconstructed due to conversions of one or both 
decay photons in the inner regions of the PHENIX detector. 
Such an effect is not included in the fastMC and is com- 
puted independently using a full simulation of the detector 
including a realistic description of the material in front of the 
EMCal. The correction factors obtained from this analysis 
are 1.067 ±0.003 for PbSc West, 1.052 ±0.004 for PbSc 
East and 1 .076 ± 0.005 for PbGl, as the material between the 
collision vertex and the EMCal is different in the east and the 
west arm and between PbGl and PbSc. 



6. Efficiency corrections of the raw r\ yields (Aii+Au) 

In the AuH-Au case, the detection efficiency is deter- 
mined using a full PISA (PHENIX Integrated Simulation 
Apphcation) GEANT3-based Monte Carlo (MC) program 
of the PHENIX detector in order to simulate the complete 
response of the calorimeter to single r| decays. The nominal 
energy resolution was adjusted in the simulation by adding 
an additional /^j-independent energy smearing of ~ 3% 
for each PbSc tower The shape, position, and width of 
the r| peak measured for all pj and centralities were thus 
well reproduced by the simulated data. The data from 
each simulated r] is embedded into real AuH-Au events and 
the efficiency for detecting the embedded r| is evaluated 
analyzing the merged events with the same analysis cuts 
used to obtain the real yields. Using this technique we 
determine efficiency corrections that account not only for the 
energy resolution and position resolution of the calorimeter, 
but also for the losses due to overlapping clusters in a real 
AuH-Au event environment. The embedding also permits 
a precise determination of the effect of edge cuts and bad 
modules. Though these effects can in principle be considered 
as geometric acceptance corrections (as done in the pH-p 
and dH-Au analyses), they depend not only on the geometry 



but also on the energy deposition of an electromagnetic 
shower in the different calorimeter towers. Lastly, in the 
full-simulation plus embedding procedure we additionally 
have control over the effects of photon conversions, as the 
GEANT simulation considers the material in front of the 
EMCal and the information whether a decay-photon converts 
is kept for evaluation in the efficiency determination. 

The input r| spectrum embedded in real events is weighted 
to match a functional form fit to the measured r| spectrum so 
that the correct folding of the r| spectrum with the resolution is 
obtained. This procedure is iterated, with the fit of the pj de- 
pendence of the input weights adjusted as the estimate of the 
efficiency correction improves, until the procedure converges 
within the nearly /77-independent statistical error of the em- 
bedded sample, approximately 3%. The final overall r| yield 
reconstruction efficiency correction factor was ~3 with a cen- 
trality dependence of <20%. The losses were dominated by 
fiducial and asymmetry cuts. 



D. r| ^ 7i"7i+7i Reconstruction 

1. Raw T| yield extraction 

The second mode of rj-meson reconstruction in PHENIX is 
via the three-body decay channel r| Tt^Ti+Ti^ with branch- 
ing ratio BR = 22.6% ± 0.4%. This mode has been used for 
the pH-p and dn-Au data, but not for the Auh-Au where the large 
detector occupancy makes the signal very difficult to extract. 
Reconstruction starts with identifying the 7t" candidates 
among the pairs of EMCal y-clusters with energy > 0.2 
GeV in the same way described in the previous section for the 
direct r| ^ yy channel. The mass of a candidate is required 
to be within two standard deviations from the peak position 
of Tt". The peak position and its width are determined by the 
71'' decay kinematics and EMCal resolution for each of the 
clusters and its position. These parameters were found to be 
consistent with the expected values. Selected candidates 
with transverse momentum pj > 10 GeV/c are assigned the 
exact mass of the meson and measured pj of the pair These 
candidates are further combined into triplets with positive 
and negative particle tracks measured by DC and PCI to have 
momentum in the range 0.2 GeV/c < pj < 4.0 GeV/c. No 
particle identification was used on the charged tracks. 

In order to extract the raw r] yields the mixed-event 
subtraction technique was not used in this case because it 
does not adequately reproduce the shape of the background 
in the real events. The most important physical reason for 
this is that there are a significant number of correlated tracks 
among the 7t+7t^ pairs coming from various heavier particle 
decays. The yield extraction was done by simultaneous fitting 
of the peak and the background in the adjacent region. The 
characteristic peak in the three -particle mass distribution 
is shown in the top panel of Fig. [T] The position of the 
peak is consistent with the nominal mass of the r|-meson 
within the statistical error of the fit shown in the figure. The 
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FIG. 7: Invariant mass distribution of pion triplets (7l*^7l+7l^) mea- 
sured in p+p collisions at ^ = 200 GeV in the range m,„,, = 0-0.9 
GeV/c^ (top panel) showing r|- and co-meson peaks. Lower panels 
showing r|-mass region m,„,, = 0.42 - 0.66 GeV/c^ demonstrate three 
different methods of the extraction of the raw yields. See text for a 
detailed explanation of each method. 



measured 8 MeV/c'^ width of the peak is narrower than in the 
r| — > YY decay channel. Unlike the YY channel where the full 
width of the peak is defined by the EMCal resolution alone, 
in Ti^Ti+Tt^ only 1/4 of the measured mass is derived from 
an EMCal-based measurement. Given the relatively low pj 
of the decay products, tracking has better resolution than the 
calorimeter. These two effects result in higher accuracy of the 
mass measurement and smaller width of the peak compared 
to the r| ^ YY analyses. 



The raw yield numbers were extracted by simultaneously 
fitting the signal with a Gaussian function and the background 
to a quadratic function. The fit was limited to the mass 
window of 510 MeV/c^ < < 640 MeV/c^. The lower 
limit is chosen to avoid the region where the 7t"7t+7i^ 
decay (branching ratio BR = 13%) yields an additional signal 
at and above m^^ = 498 MeV/c^. The upper limit is chosen 

Tz\+%- (BR = 89.1%) peak 



at 782 MeV/c^ with a width of 20-25 MeV/c^. An example 
of the fit is shown in the second panel in Fig. [T] We also 
compared the result of such a combined signal-nbackground 
fit with separate fitting of the background. For that, the region 
under the peak, 530 MeV/c^ < ms^i < 570 MeV/c^, was 
rejected from the fit and the background was approximated 
by the quadratic function. The function was interpolated and 
subtracted from the histogram in Fig. |7] (third panel). The 
histogram counts in the region initially rejected were summed 
up to calculate the yield. 

In addition, simultaneous fitting was done in the restricted 
window below <580 MeV/c^, with the background ap- 
proximated by a linear function. The same three fits were re- 
peated applying an additional condition in the analysis. Each 
charged track was required to match a hit in PC3 or in the 
EMCal in case a track missed the active area of PC3. The re- 
sultant invariant mass spectrum is shown by the lower curve in 
the top panel of Fig. [T] The amplitude of the signal is reduced 
by about a factor of two because many tracks fall outside the 
acceptance of these two detectors, but the background is also 
reduced and, more importantly, modified in its shape. The 
overall significance of the results with and without matching 
is approximately the same. Signal loss due to matching can 
be corrected with the simulation with small systematic uncer- 
tainty and the results can be compared to deduce the accuracy 
of the yield extraction procedure. Thus, for each pr point 
we obtain six statistically correlated measurements of the raw 
yields. The first measurement with its statistical error is used 
in further analysis and the variance of the six measurements 
provides the estimate of the systematic errors of the yield ex- 
traction. 



2. Acceptance and efficiency corrections of the raw r\ yields 



at a safe distance from the co 



Similar corrections as described for the YY decay channel 
need to be applied to the r| Ti^Ti+Tt^ raw yields. However, 
for the 3-pion analysis, we use the full detector simulation 
and both corrections, namely the acceptance and the efficiency 
corrections, are computed at the same time. A MC hadron de- 
cay generator was used to produce initial r|-mesons with a pr 
distribution providing satisfactory statistical significance in all 
bins after acceptance and trigger losses. The full GEANT- 
based PISA simulation was updated with the three-body de- 
cay of the r|-meson and used to decay r|-mesons. PISA also 
performs the full simulation of the PHENIX detector and gen- 
erates the response of all its subsystems up to the electronics- 
signal level, which was then processed by standard PHENIX 
reconstruction software. Special attention was paid to verify 
that the simulation code represented the real configuration of 
the detector, and that the k'^ peak parameters in the real data 
and simulation were consistent with each other The recon- 
struction of the simulated data was carried out using the same 
steps and tools as the real data. 

Figure [8] shows the combined efficiency x acceptance as a 
function of pr for the three-pion decay analysis. In order to 
compare this with the YY decay channel reconstruction effi- 
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FIG. 8: Acceptance X efficiency for the r| Tl^Tt+Jl^ as a function of 
Pj in both d+Au (dashed curve) and p+p (solid curve) collisions at 
PHENIX in Run-3. 



without applying weights. The systematic uncertainties asso- 
ciated with the measurement of the phase-space density accu- 
racy were thus obtained. This correction is shown in Fig.fTOl 
The correction factor is calculated in the range where data is 
available. 
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FIG. 9: The phase-space density of the r| 3t Jl+Jl^ decay |l65ll66l 
67] (left panel). PHENIX reconstruction and trigger efficiencies in 
p-l-p for the T] 3t"ir+7t^ decay (right panel). 



ciency one needs to multiply the acceptance curve shown in 
Fig.|6]with the obtained eT]^yy = 76% ± 3% overall efficiency 
loss. The three-body decay combined acceptance is signif- 
icantly lower than the acceptance of the yy decay channel. 
With comparable branching ratios of the two modes the re- 
sulting statistics in the three-body decay mode is expected to 
be smaller The decrease of the efficiency at high pj is due to 
the momentum cut on the to be below 4.0 GeV/c. Above 
that threshold the track sample is contaminated by products 
of in-flight decays of long-lived particles with mismeasured 
momentum. 



3. Phase-space density correction 

The r| Ti^Ti+Tt^ decay channel required an additional cor- 
rection to take into account the uneven distribution of the mo- 
menta of the three pions within the kinematically allowed re- 
gion. Such a distribution, taken from ll65il66ll67ll . is shown in 
the left panel of Fig.|9l The vertical axis is the fraction of ki- 
netic energy carried by the Ti" in the r|-meson rest frame. The 
horizontal axis shows the difference between kinetic energies 
of "KT and 71+ divided by the total in the same system. The 
left plot of Fig.|9]shows that on average the 7t"-meson caiTies 
less kinetic energy, and thus momentum, than the two charged 
7t-mesons. The right panel shows the PHENIX reconstruction 
efficiency including geometrical acceptance, high-pj- trigger 
efficiency (see next Section), and analysis cuts. The latter two 
effectively select higher momentum Tt" and lower momentum 
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FIG. 10: Phase-space density correction for p-Fp (solid curve) and 
d-l-Au (dashed curve) event samples as a function of pj . 



Trigger Corrections and Absolute Cross-Section 
Normalization 



1. Minimum-bias trigger efficiency 

The minimum bias trigger does not detect every collision; 
only a certain fraction Ztrig of the inelastic collisions and a 



7t± in the lab frame. In the Ti-meson rest frame these translate f^^^^j^^ ^ ^^^^^^ observed. The spectra 



into the effect opposite of what is shown in the left panel. In 
order to correct for that we used the following approach. The 
uniform distribution of the phase-space density produced by 
the simulated event generator was weighted according to the 
known probabilities of the 7i-meson momenta to be observed 
in the r|-meson decay. The corresponding correction was de- 
duced by comparing the reconstruction efficiencies with and 



have to be coiTected for both of these effects. The correction 
factors Ztrig/^r\, determined in t59i] for dn-Au collisions, are 
shown in Table II VI In the case of pH-p collisions, as well as in 
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TABLE IV: Correction factors (ern,;/eri), due to the efficiency of the 
MB trigger for different d+Au centralities 159,1 . 



colhsion system 


correction factor 


d+Au 0-20% central 


0.95 


d+Au 20-40% semicentral 


0.99 


d-l-Au 40-60% semiperipheral 


1.03 


d-l-Au 60-88% peripheral 


1.04 



MB dH-Au collisions, one can directly determine the inelastic 
r| cross section. Therefore, one does not apply the correction 
factors mentioned above but rather multiplies the spectra by 
the total cross section observed by the BBC, found to be 23.0 
mb ± 9.7% in Run-3 p+p collisions and 1.99 b ± 0.10 b in 
Run-3 dH-Au collisions f62!]. An additional correction has to 
be applied for the bias of the BBC to high^Bj- rj. It is found to 
be 0.79 for pH-p 1 17] and 0.94 for dn-Au 0] colUsions. 



2. High-pT Gamma-trigger efficiency 



Tt" analyses ifTil I21I1 . The Gamma triggers in PHENIX are 
implemented by adding together amplitudes in 4 x 4 adjacent 
EMCal towers during data taking and comparing them to a 
pre-set threshold. In the case of pH-p the threshold was set to 
correspond to Ey = 1.5 GeV while for dn-Au it was set at Ey 
= 3.5 GeV. In the case of AuH-Au, triggering was performed 
by a LVL2 software algorithm run over the MB-triggered 
events during data taking, such that the number of rejected 
minimum bias events were recorded. This allowed two 
different threshold triggers to be employed based on event 
centrality in AuH-Au: Ey=l.5 GeV for the 60-92% peripheral 
sample and Ey=3.5 GeV for the more central event selections. 

The trigger efficiency curves versus the energy of a single 
photon for two different threshold settings used in pH-p and 
dH-Au collisions are shown in the top panel of Fig.[TT] Based 
upon these curves, the 2-Y efficiency is calculated as for the 
previous Tt" analyses ifTil I2TI1 using the fastMC calculation. 
For this calculation, the single-photon trigger turn-on curve is 



The efficiency of the high-/??- trigger has to be studied as 
well to get r| spectra for the Gamma-triggered data at high 
transverse momenta, as previously performed for PHENIX 
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FIG. 11: Simulation result for the efficiency of the Gamma trigger 
in d-l-Au collisions. The Gamma trigger efficiency for a single y 
is shown in the top panel for p+p (solid curve) and d-l-Au (dashed 
curve). The gray band is the error of the measurement. The lower 
panel shows the recalculated trigger efficiency for YY ™d 7t*'3t+n^ 
channels for both collision systems. 
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FIG. 12: (color online) Level-2 trigger photon (top) and r| (bottom) 
efficiencies for central Au-l-Au {Ey = 3.5 GeV threshold) and periph- 
eral Au-l-Au (Ey =1.5 GeV) collisions as in Fig.[TT] For deriving the 
r| efficiencies, the histograms in the top panel were directly used, 
with the black lines denoting constant fits to the above-threshold 
asymptotic value, at >99.7% for the central trigger. The central (pe- 
ripheral) LVL2 r| sample was used only above pr > 5 (2) GeV/c. 
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represented by an integrated Gaussian for the d+Au analysis 
and by the integrated sum of two Gaussians for the p+p 
analysis. In the case of the Au+Au LVL2 triggers, the 
high-statistics measurement of the single-photon efficiency, 
which for the central trigger reaches ^100% above threshold, 
is used itself as shown in the top panel of Fig. [121 The 
derived efficiency is checked by comparing the ratio of the 
number of Tt" in MB events that carry the trigger flag to the 
number of Ti" in all MB events. In this way the normalization 
of the LVL2 data sample relative to the MB data sample is 
confirmed to be accurate to 2%. 

In the same way we determine the r| yy trigger efficiency 
which is shown in the lower panel of Figs. [TT| and [121 In 
the r| Tt^Ti+Ti^ decay channel where the statistics is very 
limited, we use the measured single-photon trigger efficiency 
curves shown in the top panel and full detector MC to deter- 
mine the efficiency of the trigger. The derived curves for pH-p 
and dH-Au are also shown in the lower panel of Fig.[TTl One 
can see that the trigger efficiency plateaus at a pr of the r| 
about twice the energy of the threshold in the case of the yy 
decay channel, but in the three-body decay mode where the 
trigger can only be fired by one of the y from n'^ yy, it re- 
quires the pt of the r] to be approximately 4 times the thresh- 
old. In central AuH-Au the r] efficiency reaches 50% (^100%) 



A final bin-shift correction is needed to take into account 
the fact that the data points of the r| spectra are plotted at the 
center of each given pr interval (bins whose width is as large 
as Apr = 2 GeV/c), which, due to the exponentially falling 
spectrum, does not represent the true physical value of the 
yield in the interval |68]. Usually, either the correction is ap- 
plied displacing the x-values horizontally (i.e. the center of 
the Pt bin is decreased) keeping their y-value, or the y-values 
are moved vertically (i.e. the yields are decreased) keeping the 
pr-values at the center of the bin. The second method (yield 
correction) is preferred here because it facilitates taking bin- 
to-bin Pt ratios of spectra (with slightly different shapes) from 
different collisions systems. The net effect of this recipe is a 
small (few %) shift downwards of the invariant r| yields in 
each Pt bin. 



F. Systematic Uncertainties 

1. T| ^ y Y analysis (p+p and d+Au) 

All systematic errors for the pH-p and the dn-Au analysis are 
summarized in Table |Vl Hereafter, the errors are categorized 
by type: 



for r| above /^j- = 5 (7 — 8 ) GeV /c, as shown in Fig. [121 bottom 
panel. The LVL2 data were used only for pT regions where 
the trigger had better than ~ 50% r] efficiency: pt > 5 GeV/c 
for the central trigger and pT > 2 GeV/c for peripheral. 



3. Cross-section normalization 

The invariant cross-sections for r| production as a function 
of pt in MB pH-p and dn-Au collisions are obtained from the 
measured number of counts in each pT bin via 

^d^a _ 1 1 1 1 N{Apt) 

d^p ~ 2npTNevt L BR Acc{pT)e{pT)ebigipT) Apt Ay ' 

(3) 

where Acc, e, and Etrfg are the acceptance, reconstruction effi- 
ciency, and trigger efficiency, respectively, determined in the 
previous Section; BR = 0.3943 ± 0.0026 is the known yy de- 
cay branching ratio of the r| meson, and L is the integrated 
luminosity obtained using the absolute inelastic cross-section 
normalization (see Section 1111 Al l. The invariant yields as a 
function of pT for a given bin in collision centrality in dn-Au 
and AuH-Au collisions are obtained via 



(4) 

I 

(A) point-to-point error uncorrelated between pT bins, 

(B) Pt correlated, all points move in the same direction but 
not by the same factor, 

(C) an overall normalization error in which all points move 
by the same factor independent of pT- 

The cross section measurement of the MB trigger has a 
type-C uncertainty of 9.7% in pH-p and 5.2% in dn-Au. All 
other systematic errors are of type B, i.e. they are pT corre- 
lated. 

The error of the raw yield (peak) extraction was estimated, 
as described in [21], calculating the error of the ratio of the 
integrals of the real and the mixed event distributions in the 
region of the background fit. The systematic error in peak 
extraction differs from the systematic error estimated for 
neutral pions in |21] because the background in the r| region 
cannot be estimated as well as the background in the Ti" 
region. This type-B error, estimated to be 4% higher than for 
pions, becomes dominant at very low transverse momenta due 
to the small S/B ratio. The error on the acceptance correction 
includes fiducial cuts on the edges of the EMCal sector as 
well as cuts around towers that have been determined to be 
hot or dead. The uncertainty in the MC (GEANT) description 
of the detector geometry is estimated varying these cuts 
slightly in the fastMC and in the embedded events (AuH-Au). 



1 ^^jVcent _ 1 1 1 jV(Ap7-,cent) 

2npT dpTdy ~ 2%ptKIT BR Acc{pT)eipT, cent) e^ngipT) AprAy 
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TABLE V: Systematic errors of the T| measurement in p+p and d+Au (Run-3) for different pj bins. The error of the peak extraction has a very 
steep slope at low pj-. 



error source 


P7-indep. 


3 GeV/c 


5 GeV/c 


10 GeV/c 


type 


peak extraction 




14.5% (p-l-p), 9.5% (d-l-Au) 


6% 


6% 


B 


geometric acceptance 




4.5% 


4.5% 


4.5% 


B 


T| reconstruction efficiency 




1.3% 


2.3% 


3.6% 


B 


global energy scale 




5.5% 


7.0% 


8.4% 


B 


energy scale linearity 




1.5% 


0.4% 


4.3% 


B 


Gamma-trigger efficiency 




9% (p-i-p), - (d-i-Au) 


0% (p-^p), 2.5% (d-HAu) 


0% 


B 


conversion correction 


2.0% 








B 


absolute cross-section normalization 


9.7% (p+p), 5.2% (d-HAu) 








C 



TABLE VI: List of systematic uncertainties in the PbSc r| measurement in Au-l-Au collisions (Run-2). Ranges generally correspond to uncer- 
tainties from the lower pj to the higher pj values of the measurement. 



error source percent error type 

raw yield (peak) extraction (point-to-point) 0-3 1 % A 

raw yield (peak) extraction (pj correlated) 10-20% B 

energy scale 3-8% B 

PID cuts 8% A 

geometric acceptance 4-2% B 

trigger efficiency 5-2% B 

reconstruction efficiency 2% A 



TABLE VII: Systematic errors of the r| 71 3t+Jl measurement in p-l-p and d-l-Au collisions (Run-3). The first number corresponds to the 
p-l-p data and the number in parentheses to d-l-Au in cases where it is different from p+p. 



error source 


PT indep. 


3 GeV/c 


4 GeV/c 


5 GeV/c 


6 GeV/c 


7 GeV/c 


8 GeV/c 


type 


EMCal geometrical acceptance 
DC -PCI acceptance 


4% (4%) 
2% (2%) 














C 
B 


acceptance variation 


0.5% (3%) 














B 


PC3-EMCal matching 


2% (2%) 














B 


TZ^ selection 




3% 


3% 


3% (3%) 


3% (3.5%) 


3% (4%) 


3% (4%) 


B 


conversion uncertainty 


3% (3%) 














C 


EMCal energy resolution 




2% 


2.5% 


3% (5%) 


4% (5%) 


5% (5%) 


5% (5%) 


B 


EMCal energy scale 




3% 


3% 


3% (4%) 


3.5% (4%) 


4% (4%) 


5% (4.5%) 


B 


Gamma-trigger efficiency 




5% 


5% 


5% (5%) 


5% (5%) 


5% (5%) 


4% (4.5%) 


B 


Gamma-trigger run-by-run variation 


4% (4%) 














B 


peak extraction in data (fit) 




10% 


13% 


20% (30%) 


13% (20%) 


23% (20%) 


30% (15%) 


A 


peak extraction in data (width) 




10% 


10% 


10% (15%) 


10% (15%) 


10% (10%) 


12% (20%) 


A 


peak extraction in simulation 


3% (5%) 














B 


branching ratio uncertainty 


1.8% (1.8%) 














C 


phase space corrections 




20% 


15% 


11% (27%) 


8% (24%) 


7% (20%) 


7% (19%) 


B 


MB trigger 


9.7% (5.2%) 














C 


trigger bias 


2.5% (1%) 














C 


total 




28% 


26% 


29% (45%) 


24% (38%) 


30% (33%) 


36% (35%) 





Those variations are found to result in differences in the 
yields of less than 5%. Different combinations of particle 
ID cuts were used in the analysis to estimate the uncertainty 
related to the photon identification. The differences among 
the various samples are less than 4% for all the different PID 
cuts for p-l-p as well as for d-i-Au reactions. The error in the 
reconstruction efficiency contains this difference. The most 
important source of uncertainty at high pr is related to the 
energy scale. The r| peak positions and widths observed in 
the data are not reproduced to better than 1.5%. An error in 
the energy scale of 1 .5% leads to an error of 4% in the yield at 
PT = 2 GeV/c and of 8% at /jj = 10 GeV/c. The eiTor of the 



high-pT trigger efficiency in p-l-p is different from in d-i-Au: it 
amounts to 7.5% at pr = 3.5 GeV/c and becomes negligible 
at = 5 GeV/c (see Section lmE2] i. 



2. T\ ^ jy analysis (Au+Au) 

The sources of systematic errors in the Au-i-Au analysis are 
listed in Table [VT] The main sources of systematic errors in 
the r| measurement are the uncertainties in the yield extraction 
(10-30%), the yield correction (^^10%), and the energy scale 
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(a maximum of '--'8%). The energy scale uncertainty is basi- 
cally the same as discussed before for the p+p and d+Au anal- 
yses. The uncertainty on the raw yield extraction was studied 
by varying the normalization region of the mixed event back- 
ground and by comparing yields extracted from 2a and 3a 
integration windows. The yields were found to vary within 
10% of the expectation for all centralities. The final results 
obtained with different PID cut combinations are found to be 
consistent within ^8%, and this was the assigned systematic 
uncertainty for the photon identification procedure. The fi- 
nal combined systematic errors on the spectra are at the level 
of ~ 10- 15% (type- A, point-to-point) and --10-15% (type-B, 
pj-correlated). 
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^7t analysis (p+p and d+Au) 



Systematic errors for the Ti'^Tt+Tt^ channel are summarized 
in Table lVlil The pH-p and dn-Au data samples have different 
systematic errors which are usually larger in dn-Au. This is 
due to the larger high-y^j- trigger threshold set during dn-Au 
data taking. The PC3-EMCal matching uncertainty is used 
to evaluate peak extraction uncertainty. The dominant sys- 
tematic uncertainties in the pH-p (dn-Au) measurement are in 
the yield extraction and the phase-space corrections, with un- 
certainties of 10-30% (10-30%) and -10% (-25%), respec- 
tively. The final combined systematic errors on the spectra are 
at the level of -30% (pH-p) and -40% (dn-Au). 

s 



IV. RESULTS AND DISCUSSION 

In this section, the fully corrected spectra for r| production 
differential in pr in pH-p, dn-Au and AuH-Au are presented, as 
well as the nuclear modification factors for dn-Au and AuH-Au 
collisions. The measured ri/Tt" ratio as a function of pj for 
the three colliding systems is presented and discussed in com- 
parison with a compilation of world data for hadron-hadron, 
hadron-nucleus, nucleus-nucleus and e+e^ collisions and to 
phenomenological (PYTHIA and "njT-scaling") expectations. 
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FIG. 13: Invariant r| cross-section as a function of transverse mo- 
mentum in p+p and d+Au collisions at y^s^^= 200 GeV measured 
in the T] ^ yy (circles) and r\ 7l^Jt+ir^ (squares) decay channels. 
The error bars of the r| — > yy are the quadratic sum of statistical and 
systematic uncertainties. The error bars (bands) of the it^Jt+ir^ spec- 
tra represent the statistical (systematic) uncertainties of the measure- 
ment. 
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A. Transverse Momentum Spectra (p+p, d+Au, Au+Au) 

The fully corrected spectra for the r| meson are shown in 
Fig.[T3]for MB events in proton-proton and deuteron-gold col- 
lisions at y/sNN = 200 GeV. The figure shows the spectra ob- 
tained in both the r| — > YY ^"d r| — > Tt^Tt+Tt^ decay channels. 
For the YY result, the error bars represent the total error, given 
by the quadratic sum of the statistical and the systematic un- 
certainties. For the pion-triplet spectra, the error bars (bands) 
represent the statistical (systematic) uncertainties. These re- 
sults agree well in spite of very different analysis approaches 
and sources of systematic uncertainties. Due to higher accep- 
tance and lower trigger threshold (see Figs. |6] [8] [TTT i. the YY 
channel has superior statistics and therefore these results alone 
are used henceforth. 



FIG. 14: Invariant r[ yields as a function of transverse momentum in 
d+Au collisions at ^/Inn = 200 GeV in four different centralities (0- 
20%, 20-40%, 40-60%, 60-88%). The error bars are the quadratic 
sum of statistical and all systematic uncertainties. For clarity, the 
data points are scaled vertically as noted in the figure. 



The invariant yields measured in four different centrality 
classes in d+Au collisions at ^/snn = 200 GeV are shown in 
Fig. [14] In Fig. [15] the fully corrected invariant spectra for 
MB and three different centrality classes in Au+Au collisions 
at y/sNN = 200 GeV are shown. The error bars represent the 
quadratic sum of the statistical and the point-to-point system- 
atic uncertainties. 
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FIG. 15: Invariant T| yields as a function of transverse momentum in 
Au+Au collisions at ^/Jjm = 200 GeV for MB and three centralities 
(0-20%, 20-60%, 60-92%). The error bars are the quadratic sum of 
statistical and point-to-point systematic uncertainties. For clarity, the 
data points are scaled vertically as noted in the figure. 
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FIG. 16: Nuclear modification factors for r| production for four d+Au 
centralities: 0-20%, 20-40%, 40-60%, 60-88%. The error bars 
(bands) around each point are the statistical (type-B systematic) un- 
certainties. The error band at R^^ = I indicates the uncertainty in 
(TdA) for esch centrality. The error box at Rj^^ = I indicates the p+p 
cross-section uncertainty of 9.7%. 



tor (Eq. |5]l. The error due to the r| reconstruction efficiency 
as well as the error due to uncertainties in the energy scale 
are very similar for the measurement of r] mesons in p+p and 
d+Au collisions as the measured data have been taken in the 
same experimental run, and they cancel almost completely in 
the ratio. 



B. Nuclear Modification Factor in d+Au, RdAu (pt) 

Medium effects in d+A collisions are quantitatively deter- 
mined using the nuclear modification factor given as the ratio 
of the measured d+A invariant yields, cP'Niia/ dprdy, over the 
measured p+p cross-sections, cficSpp/ dprdy, scaled by the nu- 
clear thickness function {Tja) in the centrality bin under con- 
sideration: 



RdAipr) 



d^NdA/dydpT 
{TdA) ■ d^Opp/dydp 



(5) 



Deviations from RdAipr) = 1 quantify the degree of departure 
of the hard d+A yields from an incoherent superposition of 
A^A^ collisions. The values of the nuclear thickness function 
for different centralities are obtained in a Glauber MC calcu- 
lation and tabulated in Table HIl The resulting RdAipr) for r| 
mesons in d+Au collisions is plotted for different centralities 
in Fig. [16] 

The data points at lower transverse momenta have large 
statistical errors. This is caused by the poor signal-to- 
background ratio of the r| peak in the sample that is not 
triggered with the Gamma trigger. The systematic uncer- 
tainties shown in the plot are computed propagating the 
experimental uncertainties in the p+p and d+Au measure- 
ments described in Section ITlIFI Some of these uncertainties 
cancel out when calculating the nuclear modification fac- 



In the case of minimum bias d+Au collisions, the nuclear 
modification factor, shown in Fig. [17] is more simply defined 
as the ratio of d+Au over p+p cross sections normalized by 
the total number of nucleons (2 • A for a d+A collision) with 
A = 197 for a gold nucleus: 



RdAipr) = 



dOdA 



2-Ada 



(6) 



pp 



All the d+Au nuclear modification factors shown in 
Figs. [16] and [it] are approximately 1 and show a very weak 
Pt and/or centrality dependence. Similar trends have been 
observed for Tt" production il7l] . As shown in the comparison 
plot of Fig. [18] the 71*^ nuclear modification factors indicate 
small shape modifications with centrality, with a possible 
Cronin enhancement on the level of 10% around 4 GeV/c 
disappearing for pj > 10 GeV/c. At high pr the MB 
result can be described well by next-to-leading-order pQCD 
calculations (xf, "69"] without implementation of the Cronin 
effect. The contribution of (anti-)shadowing effects [26, 2^] 
in the r| or 71*' production is very small, as expected for this 
kinematical region with xr =2pr/y^ ^ 0.02-0.2. 



The small role of initial-state cold nuclear effects observed 
in the mid-rapidity spectra of neutral mesons at high pr is 
also consistent with other similar observations in d+Au reac- 
tions at y/s^ = 200 GeV such as: (i) the absence of significant 
nuclear modifications in the yields of 7/^ compared to p+p 
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FIG. 17: Nuclear modification factor R^i^ for T| mesons as a function 
of pt for minimum-bias ^snn = 200 GeV d+Au collisions. The 
uncertainties are the same described in Fig. 1161 



collisions [62], and (ii) the very similar characteristics of near- 
side and away-side jet-like correlations in p+p and dn-Au iItoIi . 
Those results indicate that the nuclear medium has little in- 
fluence on the hard processes in dn-Au collisions at top RHIC 
energies and y = 0. 



C. Nuclear Modification Factor in Au+Au, Raa (pt) 

The nuclear modification factor, Raa{pt), for r| production 
in each centrality class in AuH-Au collisions is computed using 
the standard formula: 



Raa{pt) = 



(Taa) ■ d^Opp/dydp 



(7) 



where (i) the AuH-Au spectra d^N/dydpr are used in the nu- 
merator (Fig.fTsll. (ii) the pH-p invariant spectrum d^o/dydpr 
(Fig. \l3[ is used in the denominator, and (iii) (Taa) are the 
values of the average Glauber overlap function for each 
AuH-Au centrality (Table HHi. The Raa{pt) is computed taking 
the bin-to-bin ratio of AuH-Au and pH-p spectra and propa- 
gating the corresponding uncertainties. Only the acceptance 
uncertainty (^^5%) cancels in the AuH-Au/pH-p ratio of spectra. 
Figure [19] compares the nuclear modification factor for r| 
measured in central (0-20%), semicentral (20-60%) and 
peripheral (60-92%) AuH-Au collisions. The error bars are 
the total point-to-point eiTors (including type-A systematic 
and statistical uncertainties) of the AuH-Au and pH-p mea- 
surements. The error bands on the left are the uncertainties 
in (TaiiAii) for each centrality class. The error box on the 
right is the Run-3 pH-p cross-section uncertainty of 9.7%. As 
observed for high-pj- Ti^ ifTgj I21I1 . the AuH-Au r| yields are 
consistent with the expectation of independent A^A^ scatterings 
in peripheral reactions (Raa ~1) but they are increasingly 
depleted with respect to this expectation for more central 
collisions. There is no pr dependence of Raa, as seen also for 
neutral pions. 
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FIG. 18: Nuclear modification factor Rd\{pT) for ll'^ (left) and r| 
(right) production in different centrality selections and MB d+Au 
data. The bands around the data points show systematic uncertainties 
which can vary with pj (type-B errors). The shaded band around 
unity indicates the (T^ia) uncertainty and the small box on the left 
side of the data points indicates the normalization uncertainty of the 
p-l-p total inelastic cross-section. 



Figure |20] contrasts the nuclear modification factors 
measured in central AuH-Au at y/s^ = 200 GeV for r|, 
%^ lfT9l I21I1 and y ll48ll . Whereas direct photons are unsup- 
pressed compared to the scaled reference given by a NLO 
pQCD calculation IH III that reproduces the PHENIX pH-p 
Y results well ll49ll . neutral pions and r| are suppressed by a 
similar factor of ^5 compared to the corresponding cross 
sections measured in pH-p. Within the cuiTent uncertainties, 
light-quark neutral mesons at RHIC show a flat suppression in 
the range « 4 - 15 GeV/c, independent of their mass (note 
that the r| is ^4 times heavier than the n^). Those results are 
in agreement with parton energy loss calculations in a system 
with initial effective gluon densities of the order dN^ /dy ^ 
1000 (solid curve in the figure) [37]. The equal suppression 
of r| and 71*' mesons and the agreement with parton energy 
loss calculations suggest that the final fragmentation of the 
quenched parton into a leading meson occurs in the vacuum 
according to the same probabilities (fragmentation functions) 
that govern high-pj- hadroproduction in more elementary 
systems (p+p, e+e^). This conclusion is examined in more 
detail in the next two sections. 
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FIG. 19: Nuclear modification factors for r| in three Au-l-Au central- 
ities (0-20%, 20-60%, 60-92%). The errors bars are point-to-point 
uncertainties. The absolute normalization error bands at Raa = 1 are: 
(i) the uncertainties in (7^4) for each centrality (left side), and (ii) 
the p-l-p cross-section normalization uncertainty of 9.7% (right side). 
The Raa{pt) for peripheral/central Au-l-Au have been slightly dis- 
placed to the left/right (± 50 MeV/c) along the px axis to improve 
the clarity of the plot. 
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FIG. 20: (color online) R aa (p t) measured in central Au-l-Au at 
■ = 200 GeV for r|, 7l° [lallil] and for direct y |48]. The error 



bars include all point-to-point uncertainties. The error bands at Raa 
= I have the same meaning as in Fig. [19] The baseline p-l-p -^y + X 
reference used is a NLO pQCD calculation |48, 71] that reproduces 
our own data well |49], with theoretical (scale) uncertainties indi- 
cated by the dash-dotted lines around the points. The solid yellow 
curve is a parton energy loss prediction for the suppression factor 
of leading pions in a medium with initial gluon density dN^ jdy = 
1100 1371. 



D. Ratio of r] to 7i" (p-i-p, d-i-Au, Au-i-Au) 

A useful way to determine possible differences in the 
suppression pattern of and r| is to study the centrality 
dependence of the r|/7t'' ratio, Rf^iTfiij>T), in d-i-Au and 
Au-i-Au reactions and compare it with the values measured 
in more elementary systems (p-l-p, e+e^). The "world" 
\\jT^ ratio in hadronic and proton-nucleus collisions in- 
creases rapidly with pj and flattens out above pj « 2.5 



GeV/c at values R 



ri/jr" 



-^0.40 - 0.50 (see Section lIVElb . 



Likewise, in electron-positron annihilations at the Z pole 
(^ = 91.2 GeV), /?^/„o --0.5 for energetic x\ and 7i° (with 
= Phadnml Pbeam > 0.4, Consistent with the range of scaled 
momenta (z) — Phadnml Pjet considered here), as discussed in 
Section lTVE2l It is interesting to test if this ratio is modified 
in any way by initial- and/or final-state effects in d-i-Au and 
Au-i-Au collisions at RHIC energies. 

The production ratio of r| and mesons is shown in 
Fig. E] for pH-p and in Fig.|22]for dn-Au (MB and 4 centrahty 
classes). The ratio is calculated point-by-point for the d-i-Au 
measurements, propagating the corresponding errors. In the 
p-l-p case, a fit to the Ti'' spectrum Il2ll] was used. All the ratios 
are consistent with the PYTHIA ll72ll curve for p-l-p at ^/s = 
200 GeV (dashed line, see discussion in Section ITV E lb with 



an asymptotic R 



ri/jr" 



: 0.5 value. 



Figure |23] shows the Rx^/jfiipr) ratio for MB and three 
Au-i-Au centralities, obtained using the latest PHENIX Ti" 
spectra ll2ll] and removing those systematic uncertainties 
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FIG. 21: Ratio r|/ir^ measured in p-l-p collisions at y(s = 200 GeV. 
The error bars represent the point-to-point errors; the boxes represent 
the systematic uncertainties. The dashed line is the prediction of 
PYTHIA L721 for p-l-p at this cm. energy. 



which cancel in the ratio. The Ry^/Tfiipr) data for Au-i-Au is 
compared to a PYTHIA ll72ll calculation that reproduces the 
hadronic collision data well (see next Section). Within uncer- 
tainties, all the ratios are consistent with R-^/jfi ~ 0.5 (dashed 
line) and show no collision system, centrality, or pr depen- 
dence. A simple fit to a constant above pr = 2 GeV/c yields 
the following ratios: 

• R^/^o (AuH-Au cent) = 0.40 ± 0.04 (stat) ± 0.02 (syst). 
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FIG. 22: Ratio rj/Tl*' measured in different centralities in d+Au colli- 
sions at y/J^ = 200 GeV. The error bars represent all point-to-point 
uncertainties. The dashed line is the prediction of PYTHIA I72ll for 
p-l-p at this cm. energy. A few ^ri/jt°(P7') ratios have been slightly 
displaced to the left or right (± < 150 MeV/c) along the pj axis to 
improve the clarity of the plot. 



x2/ndf=0.48 

• R^/^o (AuH-Au semicent) = 0.39 ± 0.03 (stat) ± 0.02 
(syst), x^/ndf=0.26 

• Rr^/jfi (AuH-Au periph) = 0.40 ± 0.04 (stat) ± 0.02 
(syst), x^/ndf=0.42 

• R^/^ (p+p)= 0.48 ± 0.02 (stat) ± 0.02 (syst), x^/ndf = 
0.89 

• R^/^o (dH-Au)= 0.47 ± 0.02 (stat) ± 0.02 (syst), x^/ndf 
= 0.84 
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FIG. 23: Au+Au R^/jfi ratio in MB and three centrality classes (0- 
20%, 20-60%, 60-92%) as a function of pj compared to the ratio in 
d+Au and p+p collisions. The error bars include all point-to-point 
errors. The dashed line is the prediction of PYTHIA |72] for p+p at 
this cm. energy. A few R^jt^o (pt) ratios have been slightly displaced 
to the left or right (± 50 MeV/c) along the pj axis to improve the 
clarity of the plot. 



R. 



• 0.4 - 0.5 above pT-i GeV/c (Xp « 0.3). The exper- 



imental Rx^iTfiipj) ratios are also compared to PYTHIA and 
to njT-scaling expectations. PHENIX p+p, d+Au and Au+Au 
Tj/Ti'' ratios at ^Js = 200 GeV are found to be consistent with 



the obtained world data on R. 



1. r|/3t" ratio in hadron-hadron, hadron-nucleus and nuclear 
collisions (^^3- 1800 GeV) 



E. World Data on the ri/Tt" Ratio in High-Energy Particle 
Collisions 

In this last Section of the paper, we present a compilation of 
experimental r|/7t'' ratios as a function of transverse momen- 
tum, R^^^o{pt), measured in different hadronic and nuclear 
colliding systems in a wide range of center-of-mass energies 
(\A^ w 3 - 1800 GeV). The collected world data on ri/Tt" 
ratios includes: 

(i) hadron-hadron collisions (26 p+p, p+p, 7t*+p data sets), 

(ii) hadron-nucleus collisions (17 p,7t*+A sets), and 

(iii) nucleus-nucleus collisions (7 A+A data sets). 

In addition, we present also the Ry^/jfi (xp) ratio obtained from 
inclusive Ti" and r| cross-sections in e+e^ as a function of 
scaled momentum — 2 phad/V^ measured by the four LEP 
experiments at the Z pole (y/s = 91.2 GeV) . In all cases, the 
ratio R-^/jfi increases rapidly with pr (or Xp) and saturates at 



In Tables IVIIII BXl and |X] we Ust all data sets with pub- 
lished r| and Ti^ spectra and/or published R^iTfiipr) ratios 
in hadron-hadron, hadron-nucleus, and nucleus-nucleus 
collisions that we have found in the literature. Most of those 
measurements are performed around mid-rapidity. Roughly 
half of the Ry^jTioipT) listed have been directly taken from 
the original works whose references are provided in the data 
tables. A few others have been constructed by taking the 
ratio of the published Tt" and r| invariant cross-section spectra 
measured at the same In the latter case, the error in the 
ratio has been computed by adding statistical and systematic 
uncertainties quadratically. There were a few cases where the 
Pt binning of the r| spectrum did not match that of the Tt'^. 
In these cases, the Tt" spectrum was fitted with a functional 
form that reproduced the data well (usually a modified power 
law of the form discussed in ff^) and the ry/Ti^ ratio was 
then obtained by dividing the r| spectrum data points by 
the values of the Tt" function at each point. In this case, the 
error was computed by dividing the quoted r| error by the 
function value at that point. The uncertainty arising from the 
Tt" spectrum fit was obtained by computing the minimum and 
maximum ratio values at each point. Both errors were then 
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added in quadrature. 

In Tables IVIIII IIXI and |X] together with the general info 
on the collected data sets, we indicate for each measurement 
the approximate pr and xj — Ipr/y/s ranges, as well as the 
average value of R^f,fi at high pr, obtained by fitting the 
data to a constant above pj = 2 GeV/c. With the exception 
of the higher energy data iy/s^ ^ 100 GeV), most of the 
experimental ratios have been measured in a fractional 
momentum range jcj- w 0. 1 - 0.3 where the parton distribution 
functions are dominated by valence quarks (rather than 
gluons) and, hence, the produced high-/?7- Ti" and r| mesons 
come largely from q,q fragmentation. Figures |24l |25] and 
|26lshow the corresponding Ry^/jioipT) ratios for each type of 
colliding system. All the ratios show a rapid increase with pr 
and level off at R^^^o w 0.4 - 0.5 above pr ~ 3 GeV/c. No 
difference is observed for different colliding systems. The 
PHENIX p+p, d+Au and Au+Au high-/?r data presented in 
the previous section are consistent with those ratios. A fit of 
the PHENIX ratios to a constant gives R^,^o = 0.47 ± 0.03 
for both p+p and d+Au and, slightly lower but still consistent, 
^ti/jt" = 0-40 ± 0.04 for Au+Au. Together with the data points 
in Figs. I24l426l we also plot two phenomenological curves 
with PYTHIA 6.131 l>721 based on the Lund fragmentation 
model ll74ll75ll . and on mx-scaling expectations for the rj/Tt" 
ratio in p+p collisions at y/s = 200 GeV. 

a. Lund string fragmentation: The fragmentation mech- 
anism in PYTHIA is based on the phenomenological Lund 
string scheme 1, 74.. 75] which considers the color field between 
the partons to be the fragmenting entity rather than the quarks 
and gluons themselves. The string is viewed as a color flux 
tube formed by gluon self-interaction between the partons. As 
the partons move apart the potential energy stored in the string 
increases. At some point the string breaks via the production 
of new qq pairs according to the probability of a quantum- 
mechanical tunneling process, exp(— 7t j/k), which de- 
pends on the transverse-mass squared (m^ j — + pr^) and 
the string tension K w 1 GeV/fm « 0.2 GeV^. The string 
break-up process proceeds repeatedly into color singlet sys- 
tems as long as the invariant mass of the string pieces exceeds 
the on-shell mass of a hadron (each hadron corresponding to a 
small piece of string with a quark at one end and an antiquark 
at the other). At each branching, probabilistic rules are given 
for the production of flavors (uu : dd : ss = I : 1 : 0.3 in the 
default settings), spin (e.g. a 3:1 mixture between the lowest 
lying vector and pseudoscalar multiplets is used, suggested 
by spin counting arguments), and for the sharing of energy 
and momentum among the products. Regarding the latter, the 
probability that a hadron picks a fraction zof E + p^ out of the 
available E + p (p, is the momentum of the formed hadron 
along the direction of the quark q) is given by the "Lund sym- 
metric fragmentation function": 

f{z) - -z)"exp(-femTVz), (8) 

where a and b are free parameters adjusted to bring the 
fragmentation into accordance with measured LEP data, e.g. 



a = 0.3 and b = 0.58 GeV ^ are the current default values 
for PYTHIA 6.3 izl. In addition, for the flavor-diagonal 
meson states uU : dd : ss, PYTHIA also includes mixing 
into the physical mesons. This is done according to a 
parametrization, based on the mixing angles given in the 
Review of Particle Properties IstI] . In particular, the default 
choices correspond to r| = I /2{uii + dd) — 1/V2(ss) and 
ri' = 1 /2(mm + dd) + l /V2{ss). Thus, in the Tt'' - r| - r|' 
system, no account is taken of the difference in masses, an 
approximation which seems to lead to an overestimate of r| 
rates in e+e^ annihilation [103,1 . PYTHIA includes therefore 
parameters to allow an additional "tunable" suppression of r| 
andrj states. 

The PYTHIA Monte Carlo simulations of Tt" and r] 
/^T-differential cross-sections were carried out with the 
default settings. In particular, no ad-hoc suppression of 
r| was selected. Any uncertainty related to the choice of 
any (flavor-independent) settings should in principle cancel 
in the ratio of both pr spectra. As seen in Figs. 12414261 
within the (relatively large in some cases) experimental 
uncertainties, good agreement between the R-^/jfiipT) data 
and the model prediction is found for all the colliding systems 
and pt ranges, despite being at very different center-of-mass 
energies. We have also run PYTHIA at ^/s = 30 GeV as a 
reference for lower energy results, but the resulting r[/K^ 
curve, though slighdy lower at high-pj- (/J^^^oj^^) « 0.44) 
is still relatively close to the one obtained at y/s = 200 GeV. 
This is an indication that the pj dependence of the production 
mechanisms for both neutral mesons are very similar for all 
systems and cm. energies and, correspondingly, the ratio 
of /57-differential cross-sections is basically independent 
of the characteristics of the initial collision process, but 
dominated by the ratio of r| and Tt" (vacuum) FF which is 
relatively constant in this kinematic range (see discussion in 
Section lIVE2l i. 

b. mj scaling: The red shaded curve shown in Figs. |24| 
I25I a nd l26l corresponds to an empirical mx-scaling observa- 
tion IttII which assumes that the hadron differential cross- 
sections, plotted as a function of the transverse mass of the 

produced particle mj ~ \J mf + pr^, all have the same shape, 

f{mj), with an absolute normalization factor C/, which can 
vary but which is found to be the same for many species: 

E^=Ch-f{mj) (9) 



Assuming isospin symmetry for pion production, we have 
combined the measured PHENIX charged (7t++Jt^)/2 (mea- 
sured in the range pj- = 0.2 - 2.6 GeV/c) iItsIi and neutral 
(Pt =1-14 GeV/c) 1 17] pion differential nij cross-sections 
in p+p collisions and fitted them with a modified power-law 
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TABLE VIII: Hadron-hadron collisions with a published Tj/Tt" ratio and/or T| and TZ^ spectra. For each reaction we quote: the center-of-mass 
energy ^/s (and pi^h for fixed-target experiments), the pj and xt = 2pt/ \/I ranges of the measurement (the xj values are not quoted for "soft" 
spectra below pr = 1 GeV/c), and the average ti/tt" ratio above pr = 2 GeV/c obtained by fitting R-q/Tfi{pT > 2 GeV/c) to a constant. 



System 


(GeV) piah (GeV/c) pj range (GeV/c) 


xj range 


^Ti/j:»(P7->2GeV/c) 


Authors 


Collab./Exp. 


Ref. 


p+p 


1 ^ 8 

IJ.O 


1 on 16 — 94 


0.3-0.4 


0.52±0.13 


Donaldson 78 


FNAL M2 


[79] 


7t++p 




inn 1 fi ^ 


0.3-0.4 


0.49 ±0.10 


Donaldson 78 


FNAL M2 


[791 


7t^+p 


13.8 


100. 2.-3. 


0.3-0.4 


0.41 ±0.13 


Donaldson 78 


FNALM2 


[791 


7t++p 


19.4 


?nO 9—35 


0.2-0.4 


0.40 ±0.07 


Donaldson 78 


FNALM2 


[79] 


TT^+p 


1Q 4 


9nn 1^ — 4 


0.2-0.4 


0.43 ±0.04 


Donaldson 78 


FNAL M2 


179] 


p+p 


19.4 


9nn 9—35 


0.2-0.4 


0.42 ±0.04 


Donaldson 78 


FNAL M2 


[791 


p+p 




Z.OU. H". J.J 


0.2-0.4 


0.60 ±0.04 


Bonesini 88 


CERN WA70 


[801 


JI"'"+p 




Z.OU. J.J 


0.2-0.4 


0.43 ±0.05 


Bonesini 88 


CERN WA70 


[801 


7t^+p 


9'? n 


ZoU. H". J.J 


0.2-0.4 


0.57 ±0.06 


Bonesini 88 


CERN WA70 


[801 


p+p 


94 ^ 


9 S — 4 


0.2-0.3 


0.45 ±0.06 


Antille 87 


CERN UA6 


[811 


p+p 


94. 


9 S 4 


0.2-0.3 


0.48 ±0.04 


Antille 87 


CERN UA6 


[811 


p+p 


97 ^ 


400 9 — 1 fi 






Aguilar 9 1 


NA27 


[82] 


p+p 




U.O J . 


~ 0.1 -0.2 


0.55 ±0.04 


Amaldi 79 


ISR 


[831 


p+p 


30.6 


3.-4. 


0.2-0.3 


0.54 ±0.05 


Kourkoumelis 79 


ISR 


[42] 


p+p 


31.6 


530. 3.-8. 


0.2-0.5 


0.41 ±0.03 


Apanasevich 02 


FNAL 706 


[841 


p+p 


38.8 


800. 3.-8. 


0.1-0.4 


0.44 ±0.03 


Apanasevich 02 


FNAL 706 


[84] 


p+p 


52.7 


3.-6. 


0.1-0.3 


0.58 ±0.03 


Kourkoumelis 79 


ISR 


[42] 


p+p 


53. 


2.5-4. 


0.1-0.2 


0.53 ±0.03 


Akesson 85 


ISR AFS 


[851 


p+p 


53. 


2.5-4. 


0.1-0.2 


0.55 ±0.02 


Akesson 85 


ISR AFS 


[85] 


p+p 


53.2 


3.-6. 


0.1-0.2 


0.54±0.03 


Amaldi 79 


ISR 


[8J| 


p+p 


62.4 


3.-11. 


0.2-0.4 


0.55 ±0.03 


Kourkoumelis 79 


ISR AFS 


[42] 


p+p 


63. 


0.2- 1.5 




(0.07 ±0.055) 


Akesson 86 


ISR AFS 


[861 


p+p 


63. 


2.-4. 


0.06-0.13 


0.47 ±0.01 


Akesson 83 


ISR AFS 


[87] 


p+p 


200. 


2.-12. 


0.02-0.12 


0.48 ±0.03 


S.S.Adler 06 


PHENIX 


this work 


p+p 


540. 


3.-6. 


0.01 - 0.02 0.60± 0.04(5/aO ± 0. 15(^y.?f ) 


Banner 85 


CERN UA2 


[8^ 


p+p 


1800. 


12.0 


0.01 


1.02± 0. 15(^fa0 ± 0.23(i>'^f ) 


Abe 93 


CDF 


[43] 



TABLE IX: Hadron-nucleus collisions with a published Tj/Tt ratio and/or x\ and 71 spectra. For each reaction we quote: the center-of-mass 
energy ,Js^ (and pi^h for fixed-target experiments), the pj and xj = 2pt/\/s ranges of the measurement (the xj values are not quoted for 
"soft" spectra below pr = i GeV/c), and the average r|/n** ratio above pr = 2 GeV/c obtained by fitting R^/TfiipT > 2 GeV/c) to a constant. 



System 


(GeV) 


Plab (GeV/c) 


PT range (GeV/c) 


XT range 


i?^/„o(p7->2GeV/c) 


Authors 


Collab./Exp. 


Ref. 


p+Be 


19.4 


200. 


2.5-4. 


0.2-0.4 


0.28±0.15 


Povlis 83 


FNAL E629 


[89] 


p+C 


19.4 


200. 


2.-5. 


0.2-0.5 


0.58 ±0.05 


Povlis 83 


FNAL E629 


[89] 


p+Al 


19.4 


200. 


2.-3. 


0.2-0.3 


0.40±0.18 


Povlis 83 


FNAL E629 


[89] 


TT^+C 


19.4 


200. 


2.-4. 


0.2-0.5 


0.32±0.11 


Povlis 83 


FNAL E629 


[89] 


p+Be 


23.8 


300. 


2.5-5. 


0.2-0.4 


0.47 ±0.03 


Deschamps 85 


FNAL E5 15 


[901 


p+Be 


29.1 


450. 


0.1-1. 






Agakichiev 98 


TAPS/CERES 


[91] 


p+Au 


29.1 


450. 


0.1-1.2 






Agakichiev 98 


TAPS/CERES 


[21] 


p+Be 


29.1 


450. 


0.2-1.6 






Tikhomirov 95 


HELIOS 


[92] 


p+Be 


30.7 


500. 


4.-7. 


0.3-0.5 


0.40 ±0.06 


Alverson 93 


FNAL E706 


[93] 


7I"+Be 


30.7 


500. 


4.-8. 


0.2-0.5 


0.43 ±0.05 


Alverson 93 


FNAL E706 


[93] 


7I~+p 


31.1 


515. 


3.-8. 


0.1-0.5 


0.41 ±0.05 


Apanasevich 03 


FNAL 706 


[94] 


JI~+Be 


31.1 


515. 


3.-8. 


0.1-0.5 


0.48±0.01 


Apanasevich 03 


FNAL 706 


[94] 


Jt~+Cu 


31.1 


515. 


3.-8. 


0.1-0.5 


0.50±0.02 


Apanasevich 03 


FNAL 706 


[94] 


p+Be 


31.6 


530. 


3.-8. 


0.1-0.5 


0.42 ±0.01 


Apanasevich 03 


FNAL 706 


[94] 


p+Cu 


31.6 


530. 


3.-8. 


0.1-0.5 


0.42 ±0.02 


Apanasevich 03 


FNAL 706 


[94] 


p+Be 


38.8 


800. 


3.-8. 


0.1-0.4 


0.42 ±0.01 


Apanasevich 03 


FNAL 706 


[94] 


p+Cu 


38.8 


800. 


3.-8. 


0.1-0.4 


0.45 ±0.03 


Apanasevich 03 


FNAL 706 


[941 


d+Au 


200. 




2.-12. 


0.02-0.1 


0.47 ±0.03 


S.S.Adler 06 


PHENIX 


this work 
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TABLE X: Nucleus-nucleus collisions with a published T|/n" ratio and/or r| and it" spectra. For each reaction we quote: the center-of-mass 
energy ^/s^ (and pi„i, for fixed-target experiments), the px and xt = 2pjl ^JJ^ ranges of the measurement (the xj values are not quoted for 
"soft" spectra below pj =\ GeY/c), and the average r|/3t'^ ratio above pr = 2 GeV/c obtained by fitting R^^^^o {pr > 2 GeV/c) to a constant. 



System 


^A^(GeV) 


p,,h (GeV/c) 


PT rang 


e (GeV/c) 


xj range 


/?^/„o(p7->2GeV/c) 


Authors 


Collab./Exp. 


Ref. 


C-^C 


2.7 


2. 


0. 


-0.8 






Averbeck 97 


GSI TAPS 


[95] 


Ca+Ca 


2.7 


2. 


0. 


-0.7 






Averbeck 03 


GSI TAPS 


[96] 


Ni+Ni 


2.7 


1.9 


0. 


-0.7 






Averbeck 03 


GSI TAPS 


[96] 


Pb+Pb 


17.3 


158. 


0.6 


-2.6 


~ 0.1 -0.3 


0.53±0.21 


Aggarwal 00 


CERN WA98 


[97] 


S-hS 


19.4 


200. 


0.5 


-1.5 


0.1-0.2 


0.21 ±0.06 


Albrecht 95 


CERN WA80 


[98J 


S+Au 


19.4 


200. 


0.5 


-3.5 


0.1-0.3 


0.61 ±0.14 


Albrecht 95 


CERN WA80 




Au+Au 


200. 




2. 


- 10. 


0.02-0.1 


0.40 ±0.04 


S.S.Adler 06 


PHENIX 


this work 
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FIG. 24: (color online) Values of the i?T|/jt° ratios as a function of pr measured in the hadron-hadron collisions reported in Table I Villi The 
black curve is the prediction of PYTHIA |72] for the ratio in p+p at ^/s = 200 GeV, and the red shaded area indicates the empirical rax-scaling 
prescription Eq. JlOl l with fixed a= 1.2, power-law exponent n = 10. - 14., and an asymptotic RT!,^ = 0.5 ratio. 
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FIG. 25: (color online) Values of the i^ri/yt" ratios as a function of pr measured in the hadron-nucleus collisions reported in Table lIXl The 
black curve is the prediction of PYTHIA |72] for the ratio in p+p at -^J = 200 GeV, and the red shaded area indicates the empirical mj-scaling 
prescription Eq. l lIOI l with fixed a = 1.2, power-law exponent n = 10. - 14., and an asymptotic ^";^o = 0-5 ratio. 



functional form^ that reproduces the full spectra well in the 
range mt ~ 0.2 - 14 GeV/c^: 

/(mx) = {mj + ay" , with a = 1.2 and n = 10. (10) 

If one assumes that mx-scaling holds for r], then its mj = 

OT^ + pt^ spectrum can also be represented by Eqs. (|9| and 

( [Tol l (with, in principle, a different C/j) and therefore, the Tj/Tt" 
ratio as a function of pr should follow: 



mla+pT^ 



(11) 



where R'„ , o 



^^j^o — Cti/Qo is the asymptotic value of the ratio 

of Tj over Tt" for large pj-. Note that since the assumption of 
mx-scaling is that both mx -differential cross-sections have 
the same shape, the same parameters a = 1.2 and n = 10 are 



^ Note that the a and n parameters of Eq. jlOt are not independent but 
strongly correlated. They are actually related to the mean transverse mass 
of the spectrum via (mx) = 2a/{n~ 3). 



valid for both spectra as well as for the ratio (Eq. fTTT i. In all 
figures, the plotted mx-scaling curve with an asymptotic value 
of ^~y^o = 0-5 is found to be in good agreement with both the 
data and the PYTHIA predictions. We note that the agreement 
between PYTHIA and mx-scaling is not unexpected in as 
much as the Lund "fragmentation function," Eq. ([8]l, depends 
explicitly on the mj of the produced hadron. The upper red 
curve shown in all plots is that with the a and n parameters 
of Eq. ( fTTI ) that reproduce the power-law shape of the meson 
specti-a at = 200 GeV. At lower -^s, the spectra get 
increasingly steeper and a and n change accordingly (a and 
n are correlated with (mx) which itself is a logarithmically 
increasing function of \/s, i.e. (mx) =/(v^)). For illusti'ative 
purposes, we have (arbitrarily) fixed the parameter a to the 
value a = 1.2 and refitted the T\P spectra measured at different 
center-of-mass energies with n as a free exponent. With 
fixed a the corresponding values of the power-law exponent 
increase with decreasing y's as n ^ 10., 11.5, 13.5 and 14.0 
at yG = 200, 63, 27 and 13 GeV, respectively. The shaded 
red area indicates the range of expected mx-scaling ratios 
for power-law exponents n = 10.-14. The differences are 
negligible at large pj - where the r| and Ti" masses are much 
smaller than their pj and the ratio (Eq.fTTTi /?t^/„o (pr) ~ ^^/„o 
is independent of n - but become increasingly large at lower 
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FIG. 26: (color online) Values of the R-qfjfi ratios as a function of measured in the nucleus-nucleus collisions reported in Table 1x1 The 
black curve is the prediction of PYTHIA |72] for the ratio in p+p at \/s = 200 GeV, and the red shaded area indicates the empirical mx-scaling 
prescription Eq. jlOl with fixed a= 1.2, power-law exponent n = 10. - 14., and an asymptotic ^^y^o = 0-5 ratio. 



Pt (Pt ^ 3 GeVc). Furthermore, it is worth noting that in 
the low- Pt region below 1 GeV/c, the agreement between the 
data and the mx-scaling curve is not always perfect for all data 
sets, even taking into account different power-law exponents. 
This is due to the fact that at very low mj « 0.-0.4 GeV/c^, 
the pion yield rises due to contributions from multiple reso- 
nance decays and the formula (Eq.fTOli does not reproduce the 
spectral shape of the data anymore. Instead, an exponential 
behavior of the form Ed^a/d^p ~ B ■ exp{—bmT) ll99ll 
extrapolates the spectra better in the soft regime all the way 
down to mj = GeV/c^. However, for all practical purposes 
in this analysis focused on high-/??- production, we will 
consider Eq. ( fTOb (and correspondingly Eq. [TTI ) to be a good 
enough approximation. 



since our lowest pr value (pr ~ 2 GeV/c) is just in the range 
where radial flow effects start t o die out. The same holds 
true also for the recent proposal il02ll to study the r\/vP ra- 
tio as a tool to test different parton recombination scenarios 
in hadron production in nucleus-nucleus collisions. Lower-pj- 
r| measurements, which are intrinsically more difficult due to 
the reduced PHENIX acceptance and the larger yy combina- 
torial background, would be needed to better address the role 
of collective flow and/or parton recombination effects on the 
spectral shape and yields of light neutral mesons in AuH-Au 
colhsions at RHIC. 



2. Ti/ti" rafto m e+e collisions at the Z pole (yG = 91.2 GeV) 



Lastly we want to mention that in the case of nucleus- 
nucleus collisions the existence of a strong collective radial 
flow (pco// ~ 0.6 at RHIC lIlOOll ). absent in pH-p collisions, 
changes the spectral shape of different hadrons produced at 
low transverse momenta {pr < 2 GeV/c) and s hould result 
in a violation of the mx-scaling behavior llOlll . Since hy- 
drodynamical flow results in a larger boost for the (heavier) 
r| than for Ti", one expects a comparatively larger Rr^/-,^){pT) 
ratio in AuH-Au than in pH-p collisions below pj ~ 2 GeV/c. 
Unfortunately, we cannot test this assertion with RHIC data 



In this last section we are interested in determining the r|/7t'' 
ratio in an elementary colliding system such as e+e^ and com- 
paring it to the corresponding ratios obtained in hadronic and 
nuclear collisions. In e+e^ the dominant high-momentum 
hadron production mechanism is q,q fragmentation since 
gluon production (and subsequent fragmentation) occurs with 
a probability which is suppressed by a factor a^, and therefore 
plays a comparatively less significant role than in the (highest 
energy) hadronic and nuclear collisions discussed in the pre- 
vious section. Some of the experimental interest in the study 
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of r| production in e+e^ collisions was in fact triggered by 
theoretical expectations that the isoscalar mesons contained a 
significant gg component, and thus that gluon jets should ex- 
hibit an anomalously l arge tende ncy to fragment into r| and 
ri'(958) mesons 1120, [ml [T22I1 . However, this hypothesis 
was not confirmed by a de taile d analysis of the ALEPH e+e^ 
gluon fragmentation data flQ^. TablelUlists all the existmg 
measurements of inclusive 71" and r| production in e+e^ colli- 
sions at LEP at energies around the Z pole. At lower energies, 
there are several results on inclusive production in e+e^ 
but fe w r| measurements exist ( y^ = 29 and 35 GeV at SLAC 
PEP [llilllil and SLC |[nilli3,[ni[n9|], respectively), 
and we could not determine the corresponding ratios. 



TABLE XI: Experimental measurements at LEP of r\, it'^ spectra in 
e+e^ collisions at i/j' = 91.2 GeV. 



Collaboration - Year 


Particle 


Authors [Ref.] 


ALEPH 92 


Tl 


Buskulic et al. [1031 


ALEPH 96a 


n° 


Baratee?a/. [1041 


ALEPH 96b 




Barateefa/. [1051 


ALEPH 99 




Barateefa/. [1061 


ALEPH 01 


11 


Heisterefa/. [1071 


DELPHI 95 




Adam era/. [1081 


L3 91 


31° 


Adeva et al. [1091 


L3 92 


11 


Adrian! et al. [1101 


L3 94a 


7t°,ll 


Acciarrief flZ. [1111 


OPAL 98 


7tO,Tl 


Ackerstaff et al. [1121 


OPAL 00 




Abbiendi et al. [1 131 



>r|@\ls = 91.2 GeV: 
Aleph 92 
Aleph 96 
Aleph 99 
Aleph 01 
L3 92 
L3 94a 
L3 94b 
Opal 98 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 27: Differential inclusive r\ cross-section as a function of the 
scaled momentum Xp =2 p^ad / \/s measured at the Z pole by the 
three LEP experiments listed in Table IXll fitted to Eq. l |13l ) (solid 
curve). 



Figures |27] and |28] show the combined inclusive r| and Tt" 
invariant cross-sections measured as a function of the scaled 
particle momentum Xp = 2 PhadI \fs- Note that the overall r| 
and Tt" spectra have been measured in Xp ranges which are 
not completely overlapping. There are more experimental 
measurements on inclusive r| {16^) production at large (small) 

Jc^O.7 (xp <0.1). For this reason, in order to obtain the ra- 
tio of r| over 71" cross-sections, we have parametrized the T] 
cross-section as 



1 d(Sx\ 

^had dXp 



=A-{xp^br-{\-xpr 



(12) 



and taken the ratio of the individual Tt" data points over the 
resulting fit. We note that there is currently no r| FF available 
in the standard FF sets at hand in t he lit erature (BKK l,123ll . 
KKP Kretzer |!125], BFGW IHi]). Namely, the LEP 

data compiled in Fig. |27]have not been fitted and coded so 
far into any usable format that can be handled within a QCD 
collinear factorization appr oach. We are aware of only two 
works (Rolli et al. a t NLO lll27ll . and Indumathi and collab- 
orators at LO |128|]) which have tried to parametrize the r| 
FF from these data. An updated version of the r| FF would 
be useful as input to a NLO pQCD cross-section calculation 
for comparison to the results presented here and especially 
in the light of upcoming high-/??- r| asymmetry results using 
polarized beams of relevance for the proton spin program at 
RHIC [44]. Fitting all the available rj data with Eq. ((H]), we 
obtain the following empirical parametrization: 



1 t/Or| 
<^hud dxp 



: 0.0975 •(.x;p + 0.186) 



-2.953 



(1- 



^ 1.507 



, xVndf=0.37. 
(13) 



Using the fit (Eq. [T3] | and the 71" data plotted in Fig. l28]we 
have obtained the R^j.^[xp) ratio shown in Fig.|29l As seen 

for the corresponding rj/Tt*^ ratios in hadronic and nuclear col- 
lisions, at low values of (scaled) momentum the 71*^ production 
overwhelms that of r| (a significant fraction of low-energy pi- 
ons issues from decay contributions of heavier hadrons), but 



the ratio increases with Xp. 



From Xp sa 0.35 — 0.7, the ratio 



is consistent with the asymptotic ratio of 0.5 found in hadron 
and nuclear collisions (dashed curve). This Xp range corre- 
sponds to the values of fractional momenta (z) > 0.3 — 0.7 
typically carried by the leading high-pj- hadrons produced in 
high-energy h-i-p, h-i-A and A-i-A collisions |129, 130]. New 
results on inclusive r| and 71*^ production above Xp = 0.6 in 
e+e^ collisions at the B-factories (BELLE and BaBar) would 
be useful to determine whether the value of the ratio indeed 
saturates at Rj^/jfi = 0.5 or keeps increasing with Xp as sug- 
gested by Fig.|2^ 
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^7i°@\ls = 91.2 GeV: 
Aleph 96a 
Aleph 96b 
Aleph 99 
Delphi 95 
L3 91 
L3 94a 
L3 94b 
Opal 98 
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FIG. 28: Differential inclusive ir" cross-section as a function of the 
scaled momentum Xp = 2 p/^j/^/s measured at the Z pole by the four 
LEP experiments listed in Table IXll 



measured in d+Au relative to p+p (Fig.fTsTl indicates that the 
amount of nuclear shadowing and initial-state pj broadening 
is a small effect (at the 10% level) at mid-rapidity at RHIC en- 
ergies. This is in contrast with results at lower energies which 
showed a larger Cronin enhancement for high-z^j- mesons than 
observed here. One reason for the difference is likely due 
to the fact that hadron spectra at lower ,Js^ have steeper 
slopes and thus initial-state kj "kicks" produce a relatively 
larger net effect than on the harder spectra at RHIC energies. 
The unsuppressed d-i-Au yields combined with the observa- 
tion of strongly depleted yields of r| and 7i" in central Au-i-Au 
compared to binary-scaled p-i-p collisions (Fig. [T9l ) indicate 
that the suppression is a final-state effect in the hot and dense 
matter produced in the central Au-i-Au reactions. The consis- 
tent values of the rj/Tt" ratios measured at high pT in nuclear 
(Figs. |25] and |26]) as well as in more elementary p-i-p (Fig.l24b 
and e+e^ (Fig. |29l ) collisions clearly supports the idea that 
the suppression occurs at the parton level before the fragmen- 
tation of the parent quarks and gluons into a given leading 
meson. In particular, the overall agreement of the y\/Ti^ ra- 
tio measured in Au-i-Au and e+e^ collisions suggests that al- 
though the fast parent partons lose energy while traversing the 
system produced in central Au-i-Au collisions, their relative 
probability to fragment into a given meson, given by univer- 
sal fragmentation functions, is preserved as expected for final 
hadron formation in the vacuum. 




0.8 0.9 1 



FIG. 29: Ratio r|/n^ versus scaled momentum Xp = 2 pij^ij I \fs mea- 
sured in e+e^ collisions at LEP energies (Table IXH . obtained from 
the Tt" results of Fig.|28]and the r| fit, Eq. \Yi\ . The dashed line is the 
asymptotic R-^j-^ = 0.5 measured in hadronic and nuclear collisions. 



F. Summary of Experimental Results 

The studies presented here on high-/??- Tt" and r| produc- 
tion in the three colliding systems (p-i-p, d-i-Au, Au-i-Au) pro- 
vide interesting insights on initial- and final-state QCD ef- 
fects in cold nuclear matter (d-i-Au) and on the properties of 
the hot and dense medium produced in central Au-i-Au colli- 
sions. The absence of any strong deviation from the point-like 
scaling expectations for the /?7'-differential 71" and r| yields 



V. SUMMARY 

In summary, the transverse momentum spectra of r| mesons 
in the range pr = 2-12 GeV/c have been measured at mid- 
rapidity by the PHENIX experiment at RHIC in p-i-p, d-i-Au 
and Au-i-Au collisions at ^/s^ = 200 GeV. The r| mesons 
are reconstructed through their r| — > yy channel in the three 
colliding systems, as well as through the r| — > Tt^Tt+Tt^ de- 
cay mode in p-i-p and d-i-Au collisions. These data pro- 
vide additional characterization of high-pj hadroproduction 
in hadronic and nuclear collisions at RHIC energies. The 
d-i-Au yields are largely consistent with the p-i-p differential 
cross-sections scaled by the number of incoherent nucleon- 
nucleon collisions (RdA ~ 1). No pr or centrality dependence 
is observed in the nuclear modification factor within uncer- 
tainties. Such an observation indicates a null or very weak pr 
broadening and, in general, a very modest influence of cold 
nuclear matter effects, such as shadowing of parton distribu- 
tion functions, Cronin broadening, and/or hadronization by 
parton recombination, on high-/??- meson production at mid- 
rapidity at top RHIC energies. In contrast, the invariant yields 
measured in Au-i-Au are increasingly depleted with centrality 
compared to expectations from binary-scaled p-i-p collisions, 
up to a maximum factor of ^5 suppression in central colli- 
sions. The magnitude, pj and centrality dependence of the 
Au-i-Au suppression is the same for r] mesons and neutral pi- 
ons. The measured rj/jt" ratio in p-i-p, d-i-Au and Au-i-Au is 
nearly flat over pr = 2-12 GeV/c and is independent of the 
reaction centrality. A fit to a constant yields Ry^j-jfiipT) = 
0.4-0.5, in agreement with the experimental world values at 
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high pt collected here for hadron-hadron, hadron-nucleus and 
nucleus-nucleus collisions in a wide range of center-of-mass 
energies {^/sf^ 3-1800GeV), as well as at high (xp > 0.35) 
in electron-positron annihilations measured at y/s = 91.2 GeV 
at LEP. These results indicate that any initial- and/or final- 
state nuclear effects influence the production of light neutral 
mesons at large pr in the same way. The similar suppression 
pattern of r| and Ti*' mesons is consistent with the expectations 
of final-state parton energy loss in the dense medium formed 
in AuH-Au reactions. The approximately constant r|/7i'^ = 
0.40 ± 0.04 ratio measured in central AuH-Au collisions indi- 
cates that the attenuated parent partons fragment into leading 
mesons (r|,7t*') in the vacuum according to the same probabil- 
ities that govern high-pj- hadron production in more elemen- 
tary (e^e^, p+p) collisions. 



for r| production in MB pH-p and dn-Au collisions are tabu- 
lated in Table IXlJ and Table IXIIH respectively. The invariant 
dH-Au yields measured in centralities 0-20%, 20^0%, 40- 
60%, and 60-88% ai-e tabulated in Table |XlVl Finally, the 
invariant yields in AuH-Au reactions (MB and centralities 0- 
20%, 20-60%, and 60-92%) are presented in Tables|XV] The 
quoted errors are categorized by type: 

(A) is a point-to-point error uncorrected between pj bins, 

(B) is Pt correlated, all points move in the same direction 
but not by the same factor, 

(C) is a normalization error in which all points move by the 
same factor independent of pj. 

2. Nuclear modification factors (d+Au, Au+Au) 
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We report below the Raa{pt) t1 data tables for various cen- 
tralities in dH-Au and AuH-Au collisions. The errors quoted are 
the point-to-point and absolute normalization ones. Note that 
there is an additional 9.7% overall normalization uncertainty 
(Run-3 pH-p BBC error, gray box in Fig.[T9]l not tabulated. 



3. r|/n" ratios (p+p, d+Au, Au+Au) 

The ratio of r| to Ti" invariant yields in p+p, d+Au and 
Au+Au colHsion s at ,JT^ = 200 GeV at mid-rapidity are tab- 
ulated in Tables IXVIIII IXiXl and HH The data presented 
here are for minimum bias events and various centrality bins 
in d+Au and Au+Au collisions. 



APPENDIX: DATA TABLES 



1. Invariant t| cross-sections (p+p, d+Au) and yields (Au+Au) 

This Appendix collects the data tables of the pr spectra of 
r| mesons measured at mid-rapidity in p+p, d+Au and Au+Au 
collisions at ^/sTZ = 200 GeV. The invariant cross-sections 
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TABLli Xll: Inclaslic cross section measured in p+p at 200 GeV. 



PT (GeV/c) 


a (mb) 


tot. err. 


Stat. err. 


sys. err. 


error A 


error B 


error C 




2.75 


1.30e-03 


4.20e-04 


2.25e-04 


7.91e-06 





3.34e-04 


1.20e-04 


3.25 


3.78e-04 


6.81e-05 


8.09e-06 


1.37e-06 





5.79e-05 


3.48e-05 


3.75 


1.37e-04 


2.15e-05 


3.53e-06 


4.03e-07 





1.70e-05 


1.26e-05 


4.25 


5.49e-05 


8.20e-06 


1.81e-06 


1.47e-07 





6.20e-06 


5.05e-06 


4.75 


2.22e-05 


3.34e-06 


1.03e-06 


5.76e-08 





2.43e-06 


2.04e-06 


5.25 


1.08e-05 


1.70e-06 


6.82e-07 


2.83e-08 





1 20e-06 


9.90e-07 


5.75 


5.66e-06 


9.36e-07 


4.39e-07 


1.52e-08 





6.42e-07 


5.21e-07 


6.5 


2.02e-06 


3.47e-07 


1.75e-07 


5.58e-09 





2.36e-07 


1.86e-07 


7.5 


6.99e-07 


1.40e-07 


9.14e-08 


2.00e-09 





8.45e-08 


6.43e-08 


8.5 


1.81e-07 


4.92e-08 


4.04e-08 


5.35e-10 





2.26e-08 


1.67e-08 


9.5 


1.02e-07 


3.22e-08 


2.80e-08 


3.08e-10 





1.30e-08 


9.34e-09 


11 


2.21e-08 


9.24e-09 


8.52e-09 


1.35e-10 





2.94e-09 


2.03e-09 


















3.0 


7.5e-04 


2.4e-04 


1.3e-04 


2.1e-04 


l.le-04 


1.6e-04 


8.5e-05 


4.0 


8.1e-05 


2.5e-05 


1.4e-05 


2.1e-05 


1.3e-05 


1.4e-05 


9.1e-06 


5.0 


2.0e-05 


6.6e-06 


3.3e-06 


5.8e-06 


4.5e-06 


2.9e-06 


2.3e-06 


6.0 


5.8e-06 


1.7e-06 


l.le-06 


1.4e-06 


9.5e-07 


7.2e-07 


6.5e-07 


7.0 


l.Oe-06 


5.0e-07 


3.9e-07 


3.2e-07 


2.6e-07 


1.3e-07 


1.2e-07 


8.0 


4.5e-07 


2.7e-07 


2.1e-07 


1.6e-07 


1.5e-07 


5.6e-08 


5.1e-08 






TABLE Xin: Inelastic cross section measured in d+Au at , 


75;;;;;^ = 200 GeV. 






PT (GeV/c) 


a(mb) 


tot. err. 


Stat. err. 


sys. err. 


error A 


error B 


error C 




2.25 


2.06e+00 


7.49e-01 


1.30e-01 


3.37e-04 





7.30e-01 


1.08e-01 


2.75 


6.25e-01 


1.09e-01 


4.06e-02 


4.41e-05 





9.55e-02 


3.28e-02 


3.25 


1.58e-01 


2.40e-02 


1.48e-02 


7.84e-06 





1.70e-02 


8.27e-03 


3.75 


6.72e-02 


8.21e-03 


1.86e-03 


3.31e-06 





7.18e-03 


3.53e-03 


4.25 


2.57e-02 


3.18e-03 


7.56e-04 


1.28e-06 





2.78e-03 


1.35e-03 


4.75 


1.06e-02 


1.34e-03 


3.64e-04 


5.39e-07 





1.17e-03 


5.55e-04 


5.25 


4.38e-03 


5.65e-04 


1.96e-04 


2.20e-07 





4.77e-04 


2.30e-04 


5.75 


2.30e-03 


3.08e-04 


1.22e-04 


1.18e-07 





2.56e-04 


1.21e-04 


6.5 


9.20e-04 


1.27e-04 


5.11e-05 


4.86e-08 





1.05e-04 


4.83e-05 


7.5 


2.47e-04 


4.39e-05 


3.00e-05 


1.35e-08 





2.93e-05 


1.30e-05 


8.5 


9.17e-05 


1.89e-05 


1.44e-05 


5.18e-09 





1.12e-05 


4.82e-06 


9.5 


4.56e-05 


9.14e-06 


6.69e-06 


2.65e-09 





5.75e-06 


2.40e-06 


11 


1.31e-05 


2.88e-06 


2.21e-06 


7.92e-10 





1.72e-06 


6.88e-07 


















5.0 


l.le-02 


7.6e-03 


5.7e-03 


5.1e-03 


3.7e-03 


3.2e-03 


9.9e-04 


6.0 


2.7e-03 


1.6e-03 


1.2e-03 


l.Oe-03 


6.9e-04 


7.3e-04 


2.4e-04 


7.0 


5.8e-04 


4.4e-04 


3.9e-04 


1.9e-04 


1.3e-04 


1.3e-04 


5.1e-05 


8.0 


2.4e-04 


1.4e-04 


l.le-04 


8.4e-05 


6.1e-05 


5.4e-05 


2.1e-05 
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TABLE XIV: Invariant yields measured in d+Au for different centrality classes from most central (0-20%) to most peripheral (60-88%). 



Centrality 


PT (GeV/c) 


inv. yield 


tot. err. 


Stat. err. 


sys. err. 


error A 


error B 


error C 




2.25 


1.59e-03 


5.96e-04 


1.85e-04 


5.67e-04 





5.67e-04 







2.75 


5.46e-04 


1.02e-04 


5.65e-05 


8.49e-05 





8.49e-05 







3.25 


1.27e-04 


2.49e-05 


2.05e-05 


1.42e-05 





1.42e-05 







3.75 


5.84e-05 


6.89e-06 


2.43e-06 


6.45e-06 





6.45e-06 







4.25 


2.22e-05 


2.66e-06 


9.62e-07 


2.48e-06 





2.48e-06 







4.75 


9.26e-06 


1.15e-06 


4.61e-07 


1.06e-06 





1.06e-06 







5.25 


3.85e-06 


4.92e-07 


2.36e-07 


4.32e-07 





4.32e-07 





0-20% 


5.75 


1.97e-06 


2.70e-07 


1.47e-07 


2.26e-07 





2.26e-07 







6.5 


7.29e-07 


1.06e-07 


6.15e-08 


8.60e-08 





8.60e-08 







7.5 


2.26e-07 


4.57e-08 


3.64e-08 


2.76e-08 





2.76e-08 







8.5 


5.29e-08 


1.36e-08 


1.19e-08 


6.65e-09 





6.65e-09 







9.5 


3.31e-08 


8.26e-09 


7.06e-09 


4.27e-09 





4.27e-09 







11 


1.05e-08 


2.84e-09 


2.47e-09 


1.41e-09 





1.41e-09 







2.25 


1.18e-03 


4.38e-04 


1.37e-04 


4.16e-04 





4.16e-04 







2.75 


4.02e-04 


7.52e-05 


4.36e-05 


6.13e-05 





6.13e-05 







3.25 


1.02e-04 


1.93e-05 


1.59e-05 


1.09e-05 





1.09e-05 







3.75 


3.92e-05 


4.60e-06 


1.90e-06 


4.18e-06 





4.18e-06 







4.25 


1.47e-05 


1.77e-06 


7.97e-07 


1.59e-06 





1.59e-06 







4.75 


5.89e-06 


7.53e-07 


3.79e-07 


6.50e-07 





6.50e-07 







5.25 


2.99e-06 


4.01e-07 


2.35e-07 


3.25e-07 





3.25e-07 





20-40% 


5.75 


1.45e-06 


2.06e-07 


1.28e-07 


1.61e-07 





1.61e-07 







6.5 


5.44e-07 


8.24e-08 


5.41e-08 


6.22e-08 





6.22e-08 







7.5 


1.68e-07 


3.69e-08 


3.11e-08 


1.99e-08 





1.99e-08 







8.5 


5.91e-08 


1.29e-08 


1.07e-08 


7.23e-09 





7.23e-09 







9.5 


1.92e-08 


7.74e-09 


7.35e-09 


2.42e-09 





2.42e-09 







11 


1.19e-08 


2.90e-09 


2.45e-09 


1.55e-09 





1.55e-09 







2.25 


9.14e-04 


3.40e-04 


1.03e-04 


3.24e-04 





3.24e-04 







2.75 


1.85e-04 


4.31e-05 


3.25e-05 


2.83e-05 





2.83e-05 







3.25 


7.17e-05 


1.47e-05 


1.25e-05 


7.73e-06 





7.73e-06 







3.75 


2.46e-05 


3.01e-06 


1.47e-06 


2.62e-06 





2.62e-06 







4.25 


9.73e-06 


1.22e-06 


6.23e-07 


1.05e-06 





1.05e-06 







4.75 


4.19e-06 


5.67e-07 


3.27e-07 


4.63e-07 





4.63e-07 







5.25 


1.71e-06 


2.69e-07 


1.95e-07 


1.85e-07 





1.85e-07 





40-60% 


5.75 


9.73e-07 


1.71e-07 


1.33e-07 


1.08e-07 





1.08e-07 







6.5 


3.84e-07 


6.14e-08 


4.29e-08 


4.39e-08 





4.39e-08 







7.5 


1.07e-07 


3.14e-08 


2.87e-08 


1.27e-08 





1.27e-08 







8.5 


3.69e-08 


1.15e-08 


1.06e-08 


4.52e-09 





4.52e-09 







9.5 


1.41e-08 


6.31e-09 


6.06e-09 


1.77e-09 





1.77e-09 







11 


2.56e-09 


2.11e-09 


2.09e-09 


3.35e-10 





3.35e-10 







2.25 


4.07e-04 


1.54e-04 


5.17e-05 


1.45e-04 





1.45e-04 







2.75 


1.27e-04 


2.64e-05 


1.75e-05 


1.97e-05 





1.97e-05 







3.25 


3.52e-05 


7.72e-06 


6.66e-06 


3.90e-06 





3.90e-06 







3.75 


l.lle-05 


1.45e-06 


7.84e-07 


1.22e-06 





1.22e-06 







4.25 


4.55e-06 


6.09e-07 


3.40e-07 


5.06e-07 





5.06e-07 







4.75 


1.88e-06 


2.71e-07 


1.68e-07 


2.13e-07 





2.13e-07 





60-88% 


5.25 


6.34e-07 


1.21e-07 


9.76e-08 


7.10e-08 





7.10e-08 







5.75 


4.65e-07 


8.14e-08 


6.17e-08 


5.31e-08 





5.31e-08 







6.5 


2.27e-07 


4.27e-08 


3.34e-08 


2.66e-08 





2.66e-08 







7.5 


4.83e-08 


1.17e-08 


l.Ole-08 


5.87e-09 





5.87e-09 







8.5 


1.18e-08 


6.61e-09 


6.44e-09 


1.48e-09 





1.48e-09 







9.5 


6.33e-09 


3.76e-09 


3.67e-09 


8.15e-10 





8.15e-10 
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TABLE XV: Invariant yields measured in Au+Au for different centrality classes, including minimum bias (0-92%), most central (0-20%), and 
most peripheral (60-92%). 



Centrality 


Pt (GeV/c) 


inv. yield 


tot. err. 


sys. err. 


Stat. err. + error A 


error B 


error C 




2.25 


1.26e-02 


3.66e-03 


1.84e-03 


3.16e-03 


1.48e-03 







2.75 


3.90e-03 


1.12e-03 


5.78e-04 


9.61e-04 


4.66e-04 







3.25 


8.79e-04 


1.70e-04 


1.32e-04 


1.08e-04 


1.07e-04 







3.75 


2.33e-04 


5.13e-05 


3.53e-05 


3.72e-05 


2.88e-05 





0-92% 


4.50 


6.44e-05 


1.60e-05 


9.97e-06 


1.25e-05 


8.22e-06 





(MB) 


5.50 


1.14e-05 


2.57e-06 


1.82e-06 


1.81e-06 


1.52e-06 







6.50 


2.80e-06 


1.05e-06 


4.59e-07 


9.48e-07 


3.87e-07 







7.50 


9.60e-07 


2.63e-07 


1.62e-07 


2.07e-07 


1.38C-07 







8.50 


4.09e-07 


1.84e-07 


7.07e-08 


1.70e-07 


6.09e-08 







9.50 


1.51e-07 


8.25e-08 


2.67e-08 


7.81e-08 


2.32e-08 







2.25 


3.95e-02 


1.24e-02 


5.80e-03 


1.09e-02 


4.65e-03 







2.75 


1.29e-02 


3.86e-03 


1.91e-03 


3.36e-03 


1.54e-03 







3.25 


2.23C-03 


5.99e-04 


3.34e-04 


4.97e-04 


2.71e-04 







3.75 


4.53e-04 


2.07e-04 


6.88e-05 


1.95e-04 


5.61e-05 







4.50 


1.49e-04 


4.55e-05 


2.30e-05 


3.93e-05 


1.90e-05 





0-20% 


5.50 


2.74e-05 


8.33e-06 


4.37e-06 


7.09e-06 


3.65e-06 







6.50 


5.99e-06 


2.73e-06 


9.83e-07 


2.55e-06 


8.30e-07 







7.50 


2.79e-06 


6.99e-07 


4.71e-07 


5.17e-07 


4.02e-07 







8.50 


8.42e-07 


4.02e-07 


1.46e-07 


3.75e-07 


1.25e-07 







9.50 


4.18e-07 


2.18e-07 


7.40e-08 


2.05e-07 


6.43e-08 







2.25 


1.21e-02 


3.51e-03 


1.78e-03 


3.03e-03 


1.42e-03 







2.75 


2.88e-03 


8.52e-04 


4.26e-04 


7.38e-04 


3.43e-04 







3.25 


8.58e-04 


1.66e-04 


1.29e-04 


1.06e-04 


1.04e-04 







3.75 


2.54e-04 


5.40e-05 


3.85e-05 


3.78e-05 


3.14e-05 







4.50 


7.05e-05 


1.74e-05 


1.09e-05 


1.36e-05 


8.99e-06 





20-40% 


5.50 


1.57e-05 


3.36e-06 


2.51e-06 


2.24e-06 


2.09e-06 







6.50 


3.30e-06 


1.28e-06 


5.42e-07 


1.16e-06 


4.58e-07 







7.50 


1.06e-06 


2.65e-07 


1.79e-07 


1.96e-07 


1.53e-07 







8.50 


3.46e-07 


1.60e-07 


5.99e-08 


1.48e-07 


5.16e-08 







9.50 


1.54e-07 


7.97e-08 


2.72e-08 


7.49e-08 


2.36e-08 







2.25 


1.06e-03 


3.27e-04 


2.41e-04 


2.21e-04 


2.22e-04 







2.75 


3.34e-04 


1.02e-04 


7.61e-05 


6.83e-05 


7.03e-05 







3.25 


l.lle-04 


2.34e-05 


1.66e-05 


1.65e-05 


1.34e-05 





60-92% 


3.75 


4.04e-05 


9.05e-06 


6.13e-06 


6.66e-06 


5.00e-06 







4.50 


1.16e-05 


3.11e-06 


1.80e-06 


2.54e-06 


1.48e-06 







5.50 


1.67e-06 


1.07e-06 


2.65e-07 


1.04e-06 


2.22e-07 
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TABLE XVI: Nuclear modification factor Rja for r\ in d+Au collisions for different centrality classes, including minimum bias (0-88%), most 
central (0-20%), and most peripheral (60-88%). 



Centrality 


PT (GeV/c) 




tot. err. 


Stat. err. + error A 


error B 


error C 





7^ 

^. / J 


1 99 




994 


"^40 

U. JtU 


1 'KA 

U. 1 JH- 








n 9^^ 


109 


1 69 
U. 1 uz. 


117 

U. 1 1 / 




3.75 


1.24 


0.220 


0.047 


0.148 


0.137 




4.25 


1.19 


0.201 


0.052 


0.125 


0.131 




4.75 


1.21 


0.205 


0.070 


0.119 


0.133 




5.25 


1.03 


0.183 


0.080 


0.102 


0.114 


0-88% 


5.75 


1.03 


0.191 


0.097 


0.102 


0.114 


("MR 'I 


6.5 


1.16 


0.220 


0.119 


0.115 


0.128 




7.5 


0.896 


0.215 


0.160 


0.089 


0.099 




8.5 


1.28 


0.406 


0.350 


0.128 


0.142 




9.5 


1.14 


0.400 


0.356 


0.114 


0.126 




11 


1.51 


0.680 


0.635 


0.151 


0.166 




9 75 


1 15 


0.421 


931 


399 

U. JZZ, 


1 30 

U. 1 JU 




3.25 


0.922 


0.236 


1 50 

U. i JU 


0.141 


104 

U. lUt 




3.75 


1.17 


0.210 


0.057 


0.140 


0.132 




4.25 


1.11 


0.191 


0.060 


0.116 


0.126 




4.75 


1.15 


0.198 


0.078 


0.113 


0.130 




5.25 


0.982 


0.179 


0.087 


0.097 


0.111 




5 75 


056 


n 184 


103 

U. AUJ 


004 


108 

U. lUO 




6.5 


0.992 


0.199 


0.120 


0.098 


0.112 




7 5 


X87 


0.232 


0.184 


088 


101 

U. lUl 




X 5 

O.J 




n 9X5 


954 

JH- 


080 


001 

U.U7 1 




9.5 


0.894 


0.343 


0.312 


0.089 


0.101 




11 


1.30 


0.626 


0.590 


0.131 


0.148 




9 75 


1 99 


4.50 


94Q 


"^49 

U. JtZ 


1 "iO 

U. 1 J7 






I 07 


n 970 


1 fiQ 


16^ 

U. lUJ 


191 
u. IZl 




3.75 


1.14 


0.206 


0.062 


0.136 


0.129 




4.25 


1.06 


0.186 


0.067 


0.111 


0.120 




4.75 


1.05 


0.188 


0.084 


0.104 


0.120 




5.25 


1.10 


0.208 


0.112 


0.109 


0.125 


19-40% 


5.75 


1.02 


0.201 


0.120 


0.100 


0.116 




ft 5 

U.J 


1 07 


0.222 


0.141 


106 

U. lUU 


0.121 




7.5 


0.953 


0.264 


0.216 


0.095 


0.108 




8.5 


1.30 


0.425 


0.371 


0.129 


0.147 




9.5 


0.750 


0.374 


0.354 


0.075 


0.085 




11 


2.13 


0.995 


0.934 


0.214 


0.242 




7 75 


u.ou 




919 


941 


109 

U. iUZ 




^ 95 
j.z,j 


1 15 


^05 


0.202 


176 
U. 1 / u 


1 "^6 

U. 1 JU 




3.75 


1.09 


0.201 


0.071 


0.130 


0.128 




4.25 


1.08 


0.192 


0.078 


0.1 13 


0.127 




4.75 


1.14 


0.210 


0.104 


0.113 


0.135 


40-60% 

\J\J ffj 


5.25 


0.962 


0.198 


0.126 


0.095 


0.114 




5.75 


1.04 


0.234 


0.164 


0.103 


0.123 




6.5 


1.15 


0.246 


0.163 


0.114 


0.136 




7.5 


0.929 


0.314 


0.277 


0.092 


0.1 10 




8.5 


1.24 


0.490 


0.449 


0.123 


0.146 




9.5 


0.842 


0.451 


0.430 


0.084 


0.099 




11 


0.704 


0.645 


0.636 


0.071 


0.083 




9 75 


1 . JU 


515 

U.J 1 J 


"^01 

U. JUi 


"^81 

U.Jo 1 


165 

U. iUJ 




3.25 


1.30 


0.358 


0.248 


0.200 


0.158 




3.75 


1.13 


0.213 


0.085 


0.135 


0.137 




4.25 


1.16 


0.212 


0.095 


0.122 


0.141 




4.75 


1.18 


0.223 


0.119 


0.117 


0.143 


60-88% 


5.25 


0.826 


0.190 


0.138 


0.081 


0.010 




5.75 


1.15 


0.255 


0.177 


0.114 


0.140 




6.5 


1.57 


0.368 


0.269 


0.156 


0.191 




7.5 


0.969 


0.285 


0.239 


0.096 


0.117 




8.5 


0.913 


0.558 


0.538 


0.091 


0.111 




9.5 


0.874 


0.578 


0.561 


0.087 


0.106 
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TABLE XVII: Nuclear modification factor Raa{Pt) for ^1 in Au+Au collisions for different centrality classes from most central (0-20%) to 
most peripheral (60-92%). Note that there is an additional 9.7% normalization uncertainty (Run-3 p+p BBC error, gray box in Fig. [19) not 
quoted. 



(""pntralitv 


IJ I \ V.J \^ \ 1 1, f 




tot . srr. 


error C 






U. J JZ. 


997 (d'7 f\^r>\ 


C)^f\ (f\ 

U.UjU (^U.O /C } 






\J.J10 


n OQfS no 9%1 


099 (f\ 




J. / ^ 


0. 1 / o 


084 ('zLA Q%1 


019 (f\ 9^<7r>\ 








074 ni 7'^^1 


U.UIVJ v^u.o /C ) 




5.50 


0.199 


063 Gl 6%) 


014 (b 8%) 


0-20% 


6.50 


0.160 


075 (Af\ 9%1 


011 (b 8%^ 




7.50 


0.215 


067 f98 9%1 


01 5 ('6 8%") 




8.50 


0.250 


1 33 l'53 0%1 


017 ('6 8%") 




9.50 


0.222 


131 f59 3%'l 


01 5 ("6 8%^ 




2.75 


0.479 


70? f47 3%1 


037 n 8%^ 




3.25 


0.492 


0.117(23.8%) 


0.038 (7.8%) 




3.75 


0.401 


0.095 (23.6%) 


0.031 (7.8%) 




4.50 


0.446 


0.116(26.0%) 


0.035 (7.8%) 




5.50 


0.460 


0.106 (23.1%) 


0.036 (7.8%) 


20-60% 


6.50 


0.355 


0.143 (40.3%) 


0.028 (7.8%) 




7.50 


0.329 


0.095 (28.9%) 


0.026 (7.8%) 




8.50 


0.414 


0.214 (51.6%) 


0.032 (7.8%) 




9.50 


0.328 


0.193 (58.9%) 


0.026 (7.8%) 




2.75 


0.733 


0.315 (43.0%) 


0.209 (28.6%) 




3.25 


0.837 


0.211 (25.3%) 


0.239 (28.6%) 


60-92% 


3.75 


0.841 


0.208 (24.7%) 


0.240 (28.6%) 




4.50 


0.967 


0.271 (28.0%) 


0.276 (28.6%) 




5.50 


0.641 


0.415 (64.8%) 


0.183 (28.6%) 



TABLE XVIII: Ratio of r) and tC for p-l-p collisions at = 200 GeV. 



Pt (GeV/c) 


ll/Tt" 


tot. err. 


Stat. err. 


en'or A 


error B 


error C 


2.75 


0.440 


0.131 


0.076 





0.107 





3.25 


0.421 


0.041 


0.009 





0.040 





3.75 


0.446 


0.032 


0.012 





0.030 





4.25 


0.473 


0.035 


0.016 





0.031 





4.75 


0.468 


0.036 


0.022 





0.029 





5.25 


0.510 


0.045 


0.032 





0.031 





5.75 


0.561 


0.055 


0.044 





0.034 





6.5 


0.540 


0.057 


0.047 





0.033 





7.5 


0.596 


0.086 


0.078 





0.037 





8.5 


0.426 


0.099 


0.095 





0.027 





9.5 


0.588 


0.166 


0.162 





0.039 





11 


0.419 


0.164 


0.162 





0.029 
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TABLE XEX: Ratio of T| and 7i" for d+Au collisions for different centraUty classes, including minimum bias (0-88%), most central (0-20%), 
and most peripheral (60-88%). 



Centrality 


PJ \\JQ\IC) 




tot, err. 


st£it. en*. + error A 


error B 


error O 






420 


U.UJo 


028 


025 


A 
U 




L.I J 


A /l"70 
U.4/Z 


n €\AA 


A A1'^ 


A AOS 
U.UZo 


u 






Kl.DQD 


U.U4D 


A A1Q 


A ATI 


u 




J. / J 


A /m 


n (VXA 


A Al Q 
U.Uio 


A AOC 
U.UZo 


u 




'•■-ZD 


n 472 

U.'+/o 


U.UJJ 


A A17 
U.U I / 


A AOO 

U.uzy 


A 
U 




/I 7^ 




U.UjD 


A AOA 
U.UZU 


A A70 

U.uzy 


u 


V-oo/o 


J. ZD 


U.'l-DD 




A AO^ 
U.UZD 


A AOQ 
U.UZo 


u 


yNil5 ) 


^ 1^ 

J.I J 


n ^ 1 n 

U.DiU 


f\f\A'X 
U.U4J 


A AQA 


A A'31 


u 




6.5 


0.552 


0.048 


0.034 


0.033 







/.J 


478 


n n7n 
u.u/u 


064 


029 


A 
U 




O.J 


499 


n noo 
u.uyz 


087 


A A^A 
U.UJU 


A 
U 




9 5 


677 


118 


111 


041 


Q 




11 


609 


124 


119 


037 


Q 




Z.Zj 


U.:>oO 


U.UDZ 


A A/17 
U.U4/ 


A AT^ 
U.UZJ 


u 




Z. / J 




n n^;7 
u.uoz 


A r\^A 

U.UD4 


A A70 

U.uzy 


A 

u 




J. ZD 




U.UDO 


U.UOO 


022 


A 
U 




J. / J 


n 404 


A n4i 


A AOS 
U.UZo 


A A^A 
U.UOU 


A 
U 




/I 0^ 


U.J IZ 


n A/in 

U.U'+U 


U.UZD 


A All 
U.Uoi 


u 




4 7^ 


520 


A n4'3 


030 


031 


A 

u 




J. ZD 


508 


A A4Q 


037 


030 


A 

u 


U-ZU /O 


^ 7^ 

J.I J 


A <^47 


U.UDO 


045 


033 


A 

u 




6.5 


0.563 


0.063 


0.053 


0.034 







7 5 


579 


109 


104 


035 


Q 




8 5 


363 


094 


091 


022 


Q 




9.5 


0.644 


0.159 


0.154 


0.039 







11 


684 


193 


188 


041 


Q 




2.25 


0.416 


0.057 


0.05 1 


0.025 







Z. /D 


517 


A a^;q 

U.UDo 


060 


031 


A 
U 




J. ZD 


A /lO^ 
U.4ZJ 


A A77 
U.U/ / 


A A70 
U.U /z 


A at; 

U.UZD 


u 




1 7^ 
J. /D 




A A/IT 
U.U'+Z 


A A'^ 1 

U.U J 1 


A AOS 
U.UZo 


u 




't-.ZD 


A /1 7 A 


A A/1 1 


A AOQ 

u.uzy 


A AOC 
U.UZo 


u 




4 7^ 
't-. /D 


n 447 


A A4'3 


A A'3'3 
U.UJD 


A A97 
U.UZ/ 


A 
U 


on Af\c/^ 

/c 


D.ZD 


U.DZ4 


A A^7 
U.UD / 


A A/IS 


A A'31 


u 




^ 7^ 
D. /D 


n ^iQ 
u.Djy 


A A^;o 
u.uoz 


U.UDo 


A A'^O 
U.UoZ 


u 




6.5 


0.543 


0.069 


0.060 


0.033 







7 ^ 


n *;4'\ 

U.D4D 


Alio 

u. I ly 


114 


A n'i'X 


A 
U 




8 5 


546 


119 


115 


033 


Q 




9 5 


495 


206 


204 


030 


Q 




11 


1 046 


277 


269 


063 


Q 




Z.ZD 


A /ICQ 


U.UOD 


U.UDo 


A ATQ 

U.uzy 


u 




7^ 
Z. /D 


U.J /Z 




A A/;o 
u.uoy 


A A90 
U.UZZ 


A 

u 




3.25 


0.457 


0.091 


0.086 


0.027 







3.75 


0.439 


0.044 


0.035 


0.026 







't-.ZD 


U.'+DO 


A A4'3 


A A1^ 


U.UZ / 


A 

u 




A 7^ 
't-. /D 


A /I OA 


A A^l 
U.UDJ 


A A/1/1 
U.U44 


A AOQ 
U.UZV 


u 


4U-OUyO 


D.ZD 


A Af^f^ 
U.40U 


A AA/^ 
U.UDO 


A AAA 
U.UOU 


A AOC 
U.UZo 


u 




^ 7^ 

D. /D 


A <^47 


A AfiQ 
U.UOo 


082 


033 


A 
U 




6.5 


0.560 


0.078 


0.070 


0.034 







7.5 


0.528 


0.160 


0.157 


0.032 







8.5 


0.598 


0.200 


0.197 


0.036 







9.5 


0.495 


0.234 


0.232 


0.030 







11 


0.243 


0.211 


0.210 


0.015 







2.25 


0.474 


0.070 


0.064 


0.028 







2.75 


0.558 


0.090 


0.083 


0.033 







3.25 


0.500 


0.105 


0.101 


0.030 







3.75 


0.482 


0.051 


0.042 


0.029 







4.25 


0.490 


0.051 


0.042 


0.029 







4.75 


0.514 


0.062 


0.053 


0.031 





60-88 % 


5.25 


0.374 


0.068 


0.064 


0.022 







5.75 


0.569 


0.090 


0.083 


0.034 







6.5 


0.734 


0.129 


0.121 


0.044 







7.5 


0.501 


0.123 


0.119 


0.030 







8.5 


0.362 


0.216 


0.215 


0.022 







9.5 


0.469 


0.294 


0.293 


0.028 
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TABLE XX: Ratio of T] and % in Au+Au collisions for different centrality classes, including minimum bias (0-92%), most central (0-20%), 
and most peripheral (60-92%). 



Centrality 


Pt (Oev/cj 


T|/7C 


tot. err. 




Z.Zj 


Kj.jZX) 


n non /"is i 
u.uyu i^zo. 1 /c ) 




Z. /D 


A yt 1 A 


U.lzU (29.3%) 




J.ZJ 


A "XACv 

U.J4U 


U.Uol) (1 / .oyo) 






A '^AA 


r\ r\co /'^A AO/ \ 
U.Ujo (zU.Utb) 


U-y/yo 






n 'TG/ \ 

U.U83 (15. lyc) 


(JVLtJ) 






A ATI /TA ^iOf \ 

U.U/z (ZU.Dyc) 






A 1 CA 


A 1 QA /QT 1 07 \ 

U. 13U (3/.I70) 






A C/CA 


A 1 AA /IT A07 ^ 

U. iUU (i 1 » h) 




o.jU 


A /I OA 


A01A//1Q 007\ 

U.ziU (43.0%) 






A /I OA 


A KA ACS' \ 

U.ZDU ( J 1 .Uyo ) 




Z.Zj 


n /inn 

U.4UU 


A 1 0A /^n 

U. IZU ^JU.UVO ^ 




Z.I J 


A ^^A 


A 1 n(\ /QA A07 ^ 
U. i /U (3U.y%) 




1 

O.ZJ 


A Q'VA 
U.3 /U 


A 1 1 A /I A nfy^\ 
U. ilU (zy./7c) 






A O /I A 


U.ilU (45.070) 






A O/IA 


A 1 1 A /OA ^0/ \ 

U.liU (3U.0%) 






A OOA 


A 1 1 A /OO AO/ \ 

U.liU (zo.y%) 




A ^A 


A 'i^A 


U.IDU I /O) 




7.50 


0.530 


0.130 (24.5%) 




8.50 


0.470 


0.240 (51.1%) 




9.50 


0.490 


0.280 (57.1%) 




2.25 


0.360 


0.100 (27.8%) 




2.75 


0.340 


0.100 (29.4%) 




3.25 


0.370 


0.070(18.9%) 




3.75 


0.340 


0.066 (19.4%) 




4.50 


0.410 


0.096 (23.4%) 


20-60% 


5.50 


0.490 


0.096 (19.6%) 




6.50 


0.420 


0.160 (38.1%) 




7.50 


0.430 


0.098 (22.8%) 




8.50 


0.380 


0.170 (44.7%) 




9.50 


0.400 


0.200 (50.0%) 




2.25 


0.312 


0.094(30.1%) 




2.75 


0.383 


0.110(28.7%) 




3.25 


0.404 


0.081 (20.0%) 


60-92% 


3.75 


0.438 


0.093 (21.2%) 




4.50 


0.542 


0.139 (25.6%) 




5.50 


0.404 


0.257 (63.6%) 
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